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The transport coefficients of dense polymeric fluids are approximately calculated 
from the microscopic intermolecular forces. The following finite molecular weight ef- 
fects are discussed within the Polymer-Mode-Coupling theory (PMC) and compared 
to the corresponding reptation/ tube ideas: constraint release mechanism, spatial in- 
homogeneity of the entanglement constraints, and tracer polymer shape fluctuations. 
The entanglement corrections to the single polymer Rouse dynamics are shown to 
depend on molecular weight via the ratio N/Ne, where the entanglement degree of 
polymerization, Nf,, can be measured from the plateau shear modulus. Two micro- 
scopically defined non-universal parameters, an entanglement strength 1/a and a 
length scale ratio, 6 = £,p/b, where and b are the density screening and entangle- 
ment length respectively, are shown to determine the reduction of the entanglement 
effects relative to the reptation-like asymptotes of PMC theory. Large finite size 
effects are predicted for reduced degrees of polymerization up to N/Nf. < 10^. Ef- 
fective power law variations for intermediate N/N^ of the viscosity, r] ~ N^, and the 
diffusion constant, D N~y, can be explained with exponents significantly exceed- 
ing the asymptotic, reptation-like values, x > 3 and y > 2, respectively. Extensions 
of the theory to treat tracer dielectric relaxation, and polymer transport in gels and 
other amorphous systems, are also presented. 
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1. Introduction 



The transport properties of long chain polymers show characteristic phenomena attributed 
to chain "entanglements" As these effects are specific to macromolecules they appear 
connected to the internal degrees of freedom of polymers. The reptation/ tube theory has 
developed an approach to the problem of entanglements of flexible chain polymers, which 
has proved very versatile, and has found wide spread use^~^. Essential for this approach is 
the postulate of a phenomenological concept, the confining "tube" and anisotropic reptative 
motion of a polymer in the tube. 

Recently, a microscopic theoretical approach has been formulated by one of us'''^. It 
starts from approximations to the exact expressions for the microscopic forces and attempts 
to derive the dynamics of entangled polymers from the underlying equilibrium structure of 
the polymeric liquid. The connection of both theories has not been worked out compre- 
hensively, as the theoretical descriptions and some of the involved approximations strongly 
differ. However, for linear chains in 3-dimensions both theoretical approaches arrive at 
identical predictions for the exponents characterizing the asymptotic scaling of the trans- 
port properties with molecular weight. 

It is well documented that experimental studies of entangled polymeric melts and solu- 
tions find only partial agreement of predicted and measured exponents^'^. Especially, the 
long standing issue of the dependence of the shear viscosity on molecular weight should be 
mentioned; experimentally rj ~ ][£3A±o.2 observed^ , whereas both theories predict r) ~ 
asymptotically. Moreover, recent experiments on entangled polymer solutions^"^*^, and of 
the diffusion of polymer tracer chains in crosslinked gels^^'^^, find strong deviations from 
the theoretical predictions for the diffusion coefficients, D ~ M~^. Larger exponents, e.g. 
D ~ M~^'^, are reported in gels^^ and solutions, whereas in polymer melts the predicted ex- 
ponents apparently describe the limiting behavior to within experimental error At present 
the findings in gels and solutions cannot be rationalized within the reptation theory^^ even 
including non-asymptotic corrections resulting from the finite size of the tracer or matrix 
polymers^^^^^. Another set of recent intriguing experiments are the measurements of the 
dielectric relaxation times of tracer polymers in highly entangled polymeric melts^'^'^^. De- 
pending on the matrix molecular weight. Rouse, r ~ M^, reptation-like scaling, r ~ M^, 



or power law scaling with exponents in between 2 and 3 is observed. In all cases, except for 
the Rouse limit, however, the distribution of relaxation rates characterizing the dielectric 
disentanglement process is much broader than expected from the reptation results for the 
tube survival function^. 

In this paper we study whether the recent polymer mode coupling (PMC) theory^'^, 
which recovers the reptation like scalings, rj and D ~ N~'^, in the asymptotic regime, 

can account for the observed different molecular weight dependences for finite N. In partic- 
ular the PMC description of the constraint release mechanism and the constraint porosity 
corrections are analyzed within a simplified model for polymer liquids. Our goal is a unified 
understanding of all the puzzling non-asymptotic behaviors, for tracer and self diffusion 
and chain relaxation in solutions, melts and gels, within a single theoretical framework 
formulated at the molecular level. Throughout this paper extensive comparison with the 
reptation/ tube approach is drawn in order to explain the physical content of PMC theory. 
Tracer dynamical shape fluctuations are also included in our study since they are the origin 
of power law frequency behavior in the short time asymptote of the disentanglement process 
in the shear modulus or the end-to-end vector correlation function. 

In the accompanying paper^^, the theoretical description will be tested in comparisons 
with various experimental data sets. A major virtue of the PMC approach which starts from 
the microscopic force balance equations is to provide connections between the dynamics and 
the underlying equilibrium liquid and macromolecule structure. This allows independent 
theoretical or experimental information to be used in order to predict the magnitudes and 
trends of the finite size corrections and of the asymptotic prefactors. These estimates, 
and discussions of both extensions of the pure reptation/ tube picture to include finite size 
corrections and alternative non-reptation approaches to the dynamics of entangled polymers, 
will be included in the second paper. 

The contents of this paper is as follows. Section 2 presents those aspects of the theory 
which are required for the present discussion, and a careful enumeration of the necessary 
approximations. A new derivation of PMC theory based on the collective shear stress field 
as a primary slow variable is also presented. In section 3, the origins of finite size effects 
within the PMC approach are discussed, and simple models for their description are in- 



troduced. The previously obtained asymptotic predictions of PMC theory are summarized 
in section 4. A. Section 4.B develops the theoretical formulae including finite corrections 
for the diffusion constants, viscosities, and the internal and cnd-to-end-vector relaxation 
times. In section 5, model calculations for the transport coefficients, tracer and self diffusion 
constants, viscosities and dielectric relaxation times, of polymeric melts and solutions show 
the importance of finite size effects, and a discussion of their physical origins within the 
PMC description is presented. Extensions of the theory to polymer tracer diffusion through 
gels are presented and analyzed in section 6. The discussion in section 7 summarizes our 
findings and their connection to the reptation/ tube approach. Finally, we note that the 
present and following papers are rather long. We believe this is inevitable given our goals: 
(i) development of a very general theory, (ii) clear explanation of the physical and mathe- 
matical content of PMC theory and its similarities and differences with respect to reptation/ 
tube and other phenomenological approaches, and {Hi) unified, comprehensive description 
of all non-asymptotic corrections to the transport properties of melts, solutions, and gels. 

2. Theory 

A. Generalized Langevin Equations 

Use of the exact Mori-Zwanzig formalism allows one to derive a generalized Langevin equa- 
tion for the dynamics of the segments of a tracer polymer chain in a polymeric liquid^. 
In order to proceed, knowledge of the microscopic forces acting on the tracer segments is 
required. The ideality concept of Flory^^, which states that chain macromolecules in a 
dense, polymeric melt exhibit ideal, random-walk intramolecular correlations simplifies the 
intramolecular forces. The limit of large degree of polymerization, N, of the tracer, is math- 
ematically attractive since the discrete nature of the polymer segments can be neglected^'^^. 
Then the continuous Gaussian space curve description is generally assumed to be valid 
leading to F™*''^(i) = Ksd'^'Ra{t), where Ks is the entropic spring constant and Ra(^) the 
position vector of segment a on the tracer polymer. Time and segment length derivatives are 
abbreviated by 5t = ^ and d'^ = (^)^ respectively. Extensions to incorporate the discrete 
structure of the chains or non-Gaussian intramolecular correlations, like chain stiffness, have 



been studied^^'^^, but are neglected in the present work. Our approach, therefore, is hmited 
to dynamical processes acting on length scales large compared to segmental sizes. Formu- 
lating this restriction as an inequality, b ^ a is required, where b denotes the entanglement 
length (to be discussed below) and a is an effective, Gaussian-segment size determined by 
the persistence length of the semiflexible polymer^^. 

Far less is known about the time and space correlations of the intermolecular forces 
exerted by the surrounding polymeric matrix on the tracer. The rapid collective local dy- 
namical variables, including when appropriate the solvent, give rise to an instantaneous, 
uniform friction described by the monomeric friction coefficient, ^q. Yet unspecified, slow 
collective degrees of freedom lead to a non-trivial memory function Fa/sit), the autocorrela- 
tion function of the slow fluctuating intermolecular forces^ 

pN rt 

Codt Ra(i) - Ksdl Ra(t) + dp dt' r^fsit - t')dt'Rp(t') = F^(i) . (1) 

Jo Jo 

if. = ^ (2) 

r„«(«) = ^Ci(t) ■ F?(0)) , where (F«(() ■ R^(0)) = (3) 

In eqns. (1) and (3) the possible spatial anisotropy of the intermolecular forces on the seg- 
ment size length scale is neglected, and a spatially isotropic motion is assumed^'^^. Whether 
such an isotropic description is literally valid, or represents a kind of effective medium de- 
scription of cage-averaged dynamics, remains unresolved. The topology of linear chains 
determines the boundary condition of vanishing distortion at the free chain ends^^'^: 

daRait) = for a = 0, . (4) 

The initial values or equilibrium correlations again follow from the ideal Gaussian intramolec- 
ular structure^'^^: 

((R„(0)-R^(0)f)=a^|a-^|, (5) 

Special effects of semiflexibility are neglected here and eqn. (5) defines the effective segment 
size a. 



These equations have been solved and their physical predictions have been discussed by 
Schweizer and coworkers within the PMC theory based on various approximations^'^'^^"^^. 
In the present paper wc will improve upon one of the approximations, namely the frozen 
matrix assumption^, in order to study the consequences on the probe dynamics of taking 
into account more realistically the dynamics of the surrounding medium. We consider the 
tracer to be embedded in a liquid of chemically identical chain polymers of degree of 
polymerization P. The polymer segment density, Qm-i therefore is Qm = Pf/V . 

The tracer diffusion coefficient is defined in the usual manner from the motion of the 
center-of-mass, R*"'^(i) = daYiait), in the hydrodynamic, or long time, limit^'^^: 

((R°^(0 - R^^(0))2) ^ QDt for t ^ oo (6) 

Prom the previous analysis of the PMC equations it is known that they describe a dy- 
namical process which has been favorably compared to the entanglement dynamics as 
measured by diffusion experiments, rheology, pulsed field gradient NMR, and dielectric 
spectroscopy^°'^^'^^'^^'^^. The theory predicts the existence of a strongly dependent ter- 
minal or disentanglement time, , which in the present context can be defined by the final 
relaxation step in the internal, conformational dynamics of the tracer polymer^: 

^pintra^^^ • F'f'^{0)) ^ ■ Pf "(0)) c^^/'^^ for iV ^ OO . (7) 

Note that eqn. (7) holds for long times only (Markov regime). It is simplified in so far as 
the cut-off of the PMC effects below a certain entanglement length scale, 6, is not denoted 
explicitly, and also small corrections to the global modes are neglected^'^^'^^. The theory also 
predicts the existence of an entanglement plateau, Gn — QmkBT / Ng, in the shear modulus 
on intermediate times^. (The theoretical predictions for Ne are quoted below.) This result 
depends on the generally accepted, but rigorously unverified, assumption that the measured 
collective shear stress is dominated by the incoherently added, single chain contributions^'^^. 
From this amplitude of the entanglement process in the shear stress, Gn, and the final 
relaxation time of eqn. (7), follows the shear viscosity^'^^ : 

rj = GnTo + 77 - —j;^ — + V (8) 

In eqn. (8), a simple additivity assumption was used in order to describe the crossover 
from the low molecular weight, Rouse^'^^ result, 77^ = to the high molecular 



weight, asymptotic PMC result, 77 = GnTd- This assumption is expected to be rehable 
at rather high molecular weights and compares well with numerical solutions to the full 
PMC equations including crossover effects^''". 

The calculation of the tracer diffusion coefficient, and the internal final relaxation 
time, Tu , can proceed from the exact starting equations (1) to (5) only via approximations. 
Calculating the diffusion coefficient is somewhat complicated by the coupling of translational 
and rotational motions as captured in eqn. (1) in general. (See below for a further discussion 
of this point.) Following previous analysis we will neglect this coupling which arises from the 
end monomers of the tracer polymer, i.e. the boundary condition, eqn. (4), appropriate for 
hnear chains. Naturally, it is also absent for cyclic polymers. We argue that the diffusion over 
distances large compared to the size of the polymer, Ar ^ Rg — (t^Jn/Q, is determined by 
the sum of all intermolecular forces acting uniformly on the center-of-mass of the probe^^. 
Then the diffusion coefficient is determined from the memory function of the sum of all 
intermolecular forces. Reassuringly, the result for D we will get is a clear extension of 
known expressions for the diffusion coefficients of atomic or colloidal systems^^'^^ . 

The dynamics of the uniform friction exerted on the center-of-mass is unknown and 
can only be found with approximations. In the approach put forward in rcfs. 5 and 6 the 
uncontrolled but often surprisingly successful mode-coupling approximation^'^^'^^"^^ is used. 
It requires a physically motivated choice of the relevant, slow modes. We assume that the 
dynamics of the center-of-mass friction is dominated by the collective density fluctuations 
of the matrix, Cfc(i), 



In eqn. (9) the time dependent vector R!^(i) is the position vector of the segment 7 on 
the matrix polymer m. The total intermolecular fluctuating forces exerted on the tracer 
center-of-mass are therefore described by a memory function E(i) of the form^'^^: 



1 Vhr>T r d^k 

(11) 




m=l 



(9) 



and the collective monomer density of the tracer 




(10) 



where V is the volume, and the approximation results from assuming that the four point 
correlation function of the fluctuating forces is dominated by its overlap with the correlator 
of the pair variables f)_k(t)ck(t). The label Q indicates that the projected dynamics controls 
the time evolution of the fluctuating forces. The projector Q formally achieves that there is 
no further, linear coupling of the monomer coordinates to the fluctuating forces besides the 
one described explicitly^^'^^ in eqn. (1). The vertex, i.e. the normahzed overlap of the forces 
with the slow pair variables, can be calculated^ (when neglecting three point correlations. 



It is the normalized equilibrium correlation of a collective tracer variable, daF^, with the 
product of a collective matrix variable, Ck, and another collective tracer variable, dag^. 
Therefore it is natural that the total intermolecular correlation function, hk, arises. This 
aspect also holds for atomic and colloidal systems^^'^^. Whereas the collective intramolecular 
structure factor, ujk = ^^^(f?-k(0)f?k(0)), describes a single polymer, the total intermolecular 
correlation function, h{r) = hapir) — 5Q,^(r) — 1, is an inter-chain site-site radial distribution 
function correlating segments on different macromolecules. A peculiarity of polymeric tracers 
is the long-ranged spatial dependence of h{r) due to the well known correlation hole^ effect, 
h{r) ~ =ie~^"^/^9 for r :$> a. The existence of this long-ranged correlation can be argued 
on general grounds as a universal consequence of chain connectivity and interchain excluded 
volume forces. This correlation extends to the tracer size, Rg. Note that if the liquid is taken 
as a random continuum, g{r) — 1, then /ij. — > and all entanglement effects vanish. In the 
polymer reference interaction site model (PRISM) of Curro and Schweizer the correlation 
hole is rigorously recovered, but also the local structure is described reahstically^^"^^. This is 
necessary in the present context as the short-ranged excluded volume interactions dominate 
the polymer liquid structure and also determine the magnitude of the entanglement friction. 
PRISM therefore has been used to calculate the equilibrium structural correlation functions 
entering into the vertices^'^'^'^"^^. PRISM theory connects hk to the intra and collective 
structure factors, Uk and Sk respectively, where Sk = uJk + Qmhk^ and to an effective short 
ranged pseudopotential, Ck ~ —Vk/kBT, the site-site direct correlation function. 



i.e. Q = 1): 




for Q = 1 . 



(12) 



hk — CkLOkSk 



(13) 



Although the vertex is now known, the dynamics of the projected pair variables, 
Q^\i{t)c^{t) is not. The approach put forward by Schweizer^'^ uses the mode factoriza- 
tion approximation to break the average on the right hand side of eqn. (11), and interprets 
the projected single chain dynamics of the tracer as the one arising from a perturbative short 
time theory, the renormalized Rouse (RR) model (see section 2.B below) and the projected 
collective matrix dynamics as the full entangled dynamics. 

{c%{t)Q^{t)c'^Q%) ^ PiyNujf^{t)S,{t) . (14) 

Even though the intermolecular site-site correlation function, h^, is long ranged, in the 
final result for S(t) the vertex of eqn. (12) is relatively short ranged because of the anti- 
correlated normalization of the projector. 

Eqn. (15) expresses that the major contributions to the t = amplitudes of the friction 
(fluctuating forces) exerted on the tagged polymer center-of-mass arise from short-ranged, 
local intermolecular correlations. This can be seen from the wavevector dependence of the 
friction contributions in eqn. (15). Wavevectors large compared to 1/Rg, corresponding to 
local distances, are weighted most heavily as k'^uJk increases monotonically with k. How- 
ever, kRg < 1 long wavelength contributions still play a crucial role in establishing the 
total Markovian friction due to their long relaxation times. Figure ^ shows a schematic 
representation of the collective and intramolecular correlations and length scales which de- 
termine ^(t). The tracer interacts with the matrix polymers via the short ranged, effective 
potential, Ck, the direct correlation function. The spatial correlations of the fluctuating, 
intermolecular forces further propagate along the tracer and through the matrix medium. 
Intramolecular correlations along the tracer polymer are described by uj^ and are character- 
ized by two length scales, a (local) and Rg (global). The density screening length, ^p, and 
the entanglement length, b, characterize the constraining equilibrium and dynamic structure 
of the surrounding matrix. 

The diffusion coefficient, D, can now be found in the Markovian, long time limit of eqn. 

(6), 

Z,-. ^ II + E„|.^ II + (16) 



Extending previous work based on an effectively static, homogeneous matrix {Sk{t) ~ S'o), 
the present study focuses on the effects of the spatial and time dependence of the matrix 
fiuctuations, Sk{t) in cqn. (15), on the magnitudes and the molecular weight dependences 
of the diffusion constant and other transport properties. 

The calculations leading to the expressions for the internal relaxation time have 
been discussed in refs. 6, 31, and 34. In order to show the robustness of these results, 
to provide further physical interpretation, and as a basis for incorporation of all finite size 
corrections, we re-derive them arguing differently. The original aim of Kawasaki when using 
the mode-coupling approximation was to identify slowly varying, nonlinear contributions to 
the memory functions and equations of motion^^'^°'^^'^^. One experimentally observed slow 
variable in polymeric liquids is the shear stress. Its time scale, i.e. the shear disentanglement 
time, r^, grows strongly with molecular weight. It is the longest time scale, at least according 
to experimental observations, in entangled polymeric liquids not close to a demixing or order- 
disorder phase transition^^^. For example, dielectric measurements show that the internal, 
conformational relaxation time, r^, very nearly equals in melts^*^. We therefore assume 
the existence of one disentanglement time, , appearing in different collective and single 
chain properties which originates in the single chain, conformational dynamics. Thus To can 
be calculated from the disentanglement process of a (tracer) macromolecule for N = P and 
is itself affected by the matrix disentanglement process^ . 

As Td is the longest relaxation time, any variable participating in the disentanglement 
step can be considered a slow mode. Therefore, before we summarize the theoretical or 
experimental evidence about the behavior of the matrix structure factor, Sk{t), on this time 
scale, let us calculate the fluctuating force memory function, Tapit) of eqn. (3), assuming 
that it is dominated by the collective stress tensor which is known to be slow^. To that end, 
the steps in eqns. (11) to (15) will be repeated for those intermolecular friction forces which 
affect the internal conformational dynamics, replacing the collective matrix density with the 
collective matrix stress tensor. Remembering that the intermolecular forces, Fa{t), can also 
relax via the dynamics of a tracer segment /3, a new projection of Ta/s leads to: 

Here, (T^{t) is a component of the collective matrix stress tensor where one direction is 



longitudinal, i.e. parallel to the wavevector. 

a^^-Y. f UkBTli - zF-)e^''^^ . (18) 

m=l ^ 

The equilibrium averages of the components of the stress tensor are connected to the elastic 
constants or high frequency moduli of the liquid^^'^^. 

{a\ai) = ksTVGlSa, . (19) 

The formal expressions for the microscopic, elastic constants are known and connect 
to the total potential including inter- and intramolecular forces^^'^^. In order to make a 
connection to experimental measurements the elastic constants not at short, microscopic 
times but in the glass relaxation regime would be required. It is the time scale of the 
glassy or a-relaxation which is connected to the monomeric friction coefficient, Co, entering 
the Rouse model and our generahzation of it^. Except for simple liquids little theoretical 
understanding of the glassy moduli exists^^'^^. Since the exact values are not required in our 
case, we assume for simplicity: 1) The different elements of the glassy moduli are roughly 
equal, ~ G. 2) Their magnitude is G ~ QmkBT, as would be found in the Rouse 
calculation^. With these technical assumptions the vertices, which obey V^"*'"'^ = V^^'^^Sab, 
simplify to 14"'"^ = Vf. 

Non-diagonal elements in the memory function matrix, P ^ a, arise due to the connect- 
edness of the probe polymer and are a central element of the PMC approach. Monomer 
density fluctuations, g"^ = e*''^"(*\ at different tracer sites are correlated. For a Gaussian 
ideal polymer the equilibrium intramolecular correlations are well known^ and follow from 
eqn. (5): 

The standard normalization in the vertices, see eqn. (23) below, causes the appearance of 
the differential operator, (a;^^)~^, the inverse of the "matrix" eqn. (20), which is defined by: 

/ ci/3(u;f )-^u;f = 5ia - 7) , and d^u;f = ±^u;f for a = 0, TV . (21) 

Jo D 

The boundary condition simply follows from eqn. (4). The mathematical form of (c<j^^)"^ 
can easily be found, most simply from looking at the continuum limit of the known^'^^, finite 
matrix inverse cu^^. 



iu;f)-^ = 6{a-m-^A + ^)- (22) 

This leads to the following expressions for the vertices, where for simplicity the three point 
correlations are neglected, i.e. Q = 1 in eqn. (23). 

The approximations concerning the projector, Q ^ 1, and the values of the glassy moduli 
can easily be improved. Again, however, little is known about the dynamics of the four 
point correlation function {g"{t)(j^{t) ) in eqn. (17). We assume that this matrix of 
friction functions is dominated by exactly the same dynamical processes which determine 
the long time tracer diffusion coefficient. Thus, the Tapit) relaxes (in a parallel fashion) via 
probe motion described by the collective (coherent) single chain correlator evaluated in the 
renormalized Rouse model^'^, and the full collective matrix correlator, which corresponds to 
the shear modulus, Gk{t), for our present choice of slow variables. In this case in order to 
account for the matrix structure of F^^, the result which adequately captures the equilibrium, 
t — 0, correlations is: 

r cia'ci/?'(a;r')-^(a?^(t)^r(^)^T^f )K^'^)-^ - ^Gtm^iu^fru^^t) . (24) 

Jo (jJk 

Before writing down the results for the friction matrix let us note that Tai3{t) separates into 
two independent memory functions because of the structure of u!~^, eqn. (22). 

Tapit) = Co5(a - /3)(S'(t) - M{t)dl) (25) 

S'(t) is another measure of the uniform drag exerted by the matrix on all tracer segments 
and is closely related to E(t) of eqn. (15). From eqns. (17) to (25) it follows: 

E'W = (^£^ ! J^^tH^fG,(t fJ^ (26) 

3dCo J [^T^) Gk UJk 

As will become evident from section 3.B, E and E' actually agree on intermediate time 
scales. This follows as the integrand in S' differs from the one in S, eqn. (15), by a factor 
{a'^k'^ /12u!k), which equals unity for intermediate wavevectors. E and E' arc identical as long 
as the (major) contributions to the friction of the center-of-mass arise from internal modes, 
i.e. from wavevectors obeying kRg ^ 1. They differ somewhat for longer times because. 



as is evident from the different boundary conditions, the segment coordinates, Tia{t), and 
the monomer densities, Qkit), cannot be decomposed into the same normal modes^'^^ for 
finite k. This effect of the ends of the tracer chain, however, only affects the form of the 
final relaxation of the memory function. We prefer E(i) to for both simplicity reasons, 
and the fact that the well known hydrodynamic long time tail of the diffusion coefficient in 
colloidal systems is recovered^^'^^, E(t) — > for t — > oo. Very importantly, we note that 
the molecular weight scaling of the diffusion coefficient is not affected by this uncertainty 
about the final power law decay of the exact center-of -mass friction function. 

A conceptually different friction function, M[t), determines the conformational and single 
chain stress relaxations^. Eqns. (17) to (25) lead to: 

MW=<^/^(^fG.wi^. (27) 

In section 3.B the connection of eqn. (27) to the previously used expression^'^^ for M{t) is 
detailed. Let us discuss the role of the conformational memory function, M{t), in the long 
time, Markovian hmit. This is possible as the time scale of M{t) is found^ to be shorter 
than To by a factor proportional to 1/\/N. For simplicity let us also neglect the uniform 
friction contribution, i.e. S'(t), and consequently study internal, conformational dynamics 
only. Eqn. (1) then reduces to: 

Co(l - Modl)d, n^{t) - Ksdl R„(t) = F^{t) , (28) 

where 

Mo = / dt M{t) . (29) 



For large Mq, the disentanglement, or terminal relaxation time To follows from eqns. (7) 

and (28) and the use of Gaussian intramolecular forces^: 

ru = ^Mo = yMo , (30) 

where tq = a^Co/SksT connects the theoretical with the physical time. Moreover, con- 
sidering eqn. (28) for times short compared to one concludes that the conformational 
dynamics, i.e. the relaxation of the intramolecular distortions, 9^Ra(t), is arrested up to 



times short compared to Tr, . The exact non-Mar kovian analysis of refs. 6, 31, and 34 rigo- 
rizes this statement and identifies the time to be r^^ ^ tqN"^ ^ N / N^.. Solving the equation 
9^Ra(t) = Cq, the only time dependence can enter via the two constants of integration. Note, 
that they are not determined by the boundary condition, eqn. (4), as we are neglecting the 
center-of-mass motion. 

R™>^f-(t) = da'{ r da"ca" + u(t)) + v(t) . (31) 
Jo Jo 

The entanglement process in the PMC theory leads to the prediction of a very cooperative 
motion of all tracer segments. Figure |^ shows this in a discrete bead-spring model, where 
the standard discretization of 9^Ra = Ra+i — 2Rq + R^-i applies. The restriction Ha+i — 
2Rq, + Rq,-i = 0, where Cq, = is set for simplicity, constrains the conformational dynamics 
of the tracer for times shorter than the disentanglement time, ■ From the arbitrary 
(isotropic) displacements, u and v, of two randomly chosen segments the motion of all other 
monomers is determined. The PMC theory generalizes the reptation ansatz^'^, that only the 
two end segments can move freely, to a more general, molecular architecture transcending 
cooperative motion; e.g. rings and chains are expected to behave very similarly. Further 
numerical work possibly can find information about the monomer trajectories by studying 
the arrested magnitudes Cq. 

In summary, in the PMC approach two central memory functions describing dynamical 
caging of the tagged polymer by its surroundings arise in the analysis of entanglements. As- 
suming that the projected dynamics is dominated by the identical superposition of collective 
intramolecular and collective matrix structural dynamics, both memory functions are closely 
connected. For linear chains, in contrast to cyclic polymers, a technical complication arises 
from the end monomers but will be neglected as argued above and elsewhere^'^^. The close 
connection of T,{t) and M{t) immediately leads to the prediction that from eqns. (15) and 
(27), and (reasonable) results for the collective matrix dynamics (see below), the asymptotic 
Stokes-Einstein ratio obeys^'^^'^^: 

rjD/Rl ~ for ^ oo . (32) 

This PMC prediction agrees with the reptation result and shows that the dynamics of 
entangled polymers differs strongly from the Rouse modeP'^, rj^D^ ~ N^. In this paper we 



will analyze finite- A?" corrections to D, rj, and in order to study whether experimentally 
observed deviations from asymptotic scahng laws are described by the finite size static and/ 
or dynamic corrections of the PMC theory. 

B. Renormalized Rouse Model 

The collective, single tracer projected dynamics entering the friction functions, E(i) and 
M(t) eqns. (15) and (27), is required. Corrections to the Rouse model on short time scales 
have been worked out and discussed in references 6 and 32. Wc will use the Markovian results 
of this Renormalized Rouse (RR) theory as we are interested in the transport properties of 
the PMC model only. In the large N limit (effectively frozen matrix) , the decay of the PMC 
friction is dominated by the collective tracer dynamic structure factor evaluated in the RR 
model^'^^: 

c.rW=-^.exp[-^|^], (33) 

where is the friction coefficient of the RR model. It results from a dynamically per- 
turbativc calculation of the intermolecular friction and, for large N, exceeds the monomeric 
friction coefficient^^ Cq: 

C^^ = Co{S/27)g^d'gl / dkeculSk ^ U— for TV ^ oo , (34) 

where, for simplicity we drop the regular, molecular weight independent contribution (Co) to 
^RR g^j^^ Yodk at times long compared to the Rouse time, = Tq{N/t^Y, where attains its 
Markovian value^'^^. gd is the value of the intermolecular segment-segment pair correlation 
function, g{r), at the excluded volume diameter d, i.e. at the distance of closest approach or 
contact. The prefactor of the asymptotic limit, (^^^ ~ ^/N ■, defines the entanglement degree 
of polymerization, N^. in the theory. It is defined to be the N where in a perturbative, 
crossover calculation the increase of the friction due to the slow tracer-matrix interactions 
equals the instantaneous friction modeled by Co in the Rouse picture^^'^^. A^e defines the 
crossover length scale, h = ay/N^, and the Rouse mode index, pe = which separates 
local. Rouse dynamics from the strongly entanglement affected, more collective modes. In 
the PMC theory the entanglement effects on local length scales are explicitly removed^'^^. 
As a cut-off length or Rouse mode index the quantities b or pe are chosen, respectively. 



For later reference we quote the result for Ne following from eqn. (34) when 

Setting29'31,32,34 ^RR(^ ^ ^ 

TVe = [V37r(16/9)^„(73r-6^2^o]-' , (35) 

where F = cr/d is an effective chain aspect ratio^^'^^, set equal to unity for simplicity through- 
out this and the following paper. 

Without repeating the discussion of the RR model from references 6 and 32, we recall 
that the build-up of friction is evaluated for times short compared to the dynamics of 
the matrix. Thus, the assumption in the RR model of a frozen matrix, Sk{t) — Sk, is 
appropriate for entangled polymeric liquids only. Extensions of the theory to the tracer 
dynamics in short, unentangled polymer solutions and melts require a different treatment 
than presented here. The RR results, eqns. (33) and (34), can be viewed as a weak coupling 
and short time limit to the E-memory function of the PMC description. Namely, the effective 
potential or direct correlation function is estimated from the excluded volume interactions, 
c(r) = —V{r)/kBT = —{4:7i/3)d^gdS{r), and the collective dynamics of the tracer is taken 
from the Rouse model and is assumed to be rate determining for the fluctuating force 
relaxation. 

3. Finite Size Eflfects in PMC Theory 
A. Tracer Shape Fluctuations 

The relaxation of the entanglement constraints and their resulting memory functions, E(i) 
and M(t), progresses via parallel dynamical processes or channels. If any of these channels 
were ineffective, i.e. its contribution to the friction did not relax, the entanglements could 
only relax via the other channels and would be slowed down. Besides the matrix relaxation, 
which will be discussed in the next section, the tracer dynamics as described by eqn. (33) 
opens two channels for fluctuating force relaxation. A more collective one at small wavevec- 
tors is associated with probe center-of-mass translation, and a more local decay channel at 
larger wavevectors which can be identified with probe shape fluctuations. The first process 
relaxes the fluctuating forces via the coherent motion of all tracer segments, i.e. the center- 
of-mass dynamics of the RR model, ujf^{t) ~ ]S[e~(''BTk'^t/NC^^) ^ ^ Presumably 



this decay channel is the PMC (isotropic) analog of a coherent (anisotropic) reptative mo- 
tion. The local probe shape fluctuations for kRg 3> 1 arise from the continuous spectrum 
of internal modes. Eqn. (33) approximates those with the well known first cumulant ex- 
pression advocated by Akcasu and coworkers^^ which correctly describes the time-spatial 
correlations, i.e. the scaling with kH in the case of the Rouse modeP'^. 

The shape, or internal mode, fluctuations lead to a speeding up of the dynamical re- 
laxation of fluctuating (entanglement) forces associated with dynamically correlated pro- 
cesses on length scales <^ Rg. It is instructive to neglect the shape fluctuation mechanism 
and extend the diffusive center-of-mass dynamics, <l>^(i, /c) = ^-{kBTkH/NQ^'^) ^ arbitrary 
wavevectors and compare the predicted behavior with the full PMC and the reptation/ tube 
results. In order not to vary the wavevector dependent contributions to the entanglement 
friction functions this change in the decay rates has to be accompanied by an appropriate 
change in the entanglement amplitudes. Adopting the frozen matrix approximation for sim- 
plicity these adjustments lead, starting from eqns. (15) and (27), to the following memory 
functions: 



By construction this change does not lead to changes in the asymptotic scaling of the trans- 
port coefficients, but has interesting effects for the short time asymptote of the disentan- 
glement process. Note also that the identical A^-scaling for the internal dynamics, i.e. 
Mq ~ Mq ~ N^, would be obtained even without adjusting the entanglement amplitudes 
but only suppressing the shape fiuctuation contribution to force relaxation. It is only in 
the center-of-mass friction function, E(i) eqn. (15), that neglecting the shape fluctuations 
without correcting the amplitudes would lead to a non-reptation like Markovian result. 



Let us recall from the reptation/ tube theory that the polymer end-to-end vector cor- 
relation function, (P(t) ■ P(0)), and the shear modulus, G(t), are proportional to the same 
function '0(t), the tube survival function^. Therefore, they exhibit the following asymptote 
in the reptation theory: 




(36) 




(37) 



G"{u)/Gn oc (P"(a;) • P(0))/i?J oc {utqN^)-^''^ for utqN^ > 1 reptation . (38) 

In PMC theory the end-to-end vector correlation function and the shear modulus are not 
rigorously proportional to each other in general. Elsewhere^°'^^'^^, it was found that in the 
frequency window -C a; <^ 1/{tqN^) the following connections hold: 

M"(uj) 

G'{uj) ~ On , and G"{uj) ~ On ^ ' , (39) 

(P'(a;) • P(0)) ~ Rl , and (P"(a;) • P(0)) ~ i^^^^^^j^f^ , (40) 

where the single/ double primes denote storage/ loss functions. The simplified diffusive 
center-of-mass relaxation of the entanglement constraints as described by the memory func- 
tions E and M, eqns. (36) and (37), leads to the identical results as the reptation tube model, 
i.e. eqn. (38). In the frequency window 1/r^^ <^ a; <C 1/t^, considering the reptative-like 
coherent decay channel only, then eqns. (34,36,37,39,40) immediately lead to eqn. (38). 
These reptation-like results in intermediate frequency windows have to be contrasted with 
the correct PMC results^'^'^^'^^ which, in the specified frequency window and using the frozen 
matrix approximation, follow from eqns. (15,27,33,39,40): 

G"{uj) ~ Gn{ujtoN^/Y^/^ , (41) 

(P"(a;) ■ P(0)) ~ i?^(a;roiV23/6)-3/8 . (42) 

Without derivation let us recall the results for higher frequencies but still below the crossover 
frequency to the Rouse behavior^O'^^'^^: G"{uj) - G'iv(a;roiV^6/3)-3/i6 (P"{uj) ■ P(0)) - 

i?2(curoAr40/9)-9/32. 

Reptation and PMC theory therefore describe similar friction contributions acting on the 
tracer polymer. The asymptotic A^-scaling of the transport coefficients agrees, but a different 
mode spectrum or intermediate time dependence arises from the treatment of the internal 
shape fluctuations in the PMC approach. Neglecting the shape fluctuations PMC theory 
recovers the reptation results as both theories then assume that entanglements only relax 
via the coherent center-of-mass diffusion of the tracer polymer. Tracer shape fluctuations 
in the PMC equations do not lead to changes in the asymptotic, Markovian results for the 



conformational dynamics, but do affect the initial decay in the final disentanglement process. 
It is just in this intermediate frequency window where the exponents and the overall shapes of 
the shear moduli predicted by PMC theory arc in much better agreement with experimental 
measurements than the reptation/ tube predictions^^'^'^'^^^^^. Note that shallow slopes, 
G"{ijj) ~ are found in the shear moduh of various polymeric systems^°~^^, whereas 

in dielectric loss^°'^^~^^ often somewhat higher exponents are measured, e"{u) ~ to 
~ a; -0-33. 

The extensions of reptation including contour fluctuations lead to shallower slopes in 
G"{ijj) than pure reptation and also model the anomalous molecular weight dependence of 
the shear viscosity^^'^^. However, in our opinion these approaches do not accurately reproduce 
the observed power law N and uj scalings. Moreover, recent dielectric tracer experiments of 
Adachi and coworkers^^'^^ find on the one hand a reptation-like scaling of the final relaxation 
time. To ~ M^, for large polymer tracers in entangled polymeric matrices, but in the identical 
systems a shallow loss spectrum, e"{ijj) ~ to ~ uj~^-^^. They conclude this disproves 

the idea that the anomalous exponents in the viscosity versus molecular weight scaling and 
in the dielectric loss spectrum are connected. Such a conclusion is in good agreement with 
PMC theory, where two very different physical effects arc the source of these behaviors. 
As discussed later in this section, within PMC theory tracer shape fluctuations lead to 
anomalous frequency power laws for the initial stages of the disentanglement process, and, 
finite size effects of the matrix constraints lead to non-reptative scaling of the transport 
coefficients. 

B. Matrix Constraint Porosity and Constraint Release 

Contributions to the memory functions E(t) and M(t), eqns. (15) and (27), are characterized 
by spatially varying or wavevector dependent amplitudes and characteristic times. Within 
the PMC description one may view the static contributions or amplitudes in the mode 
coupling vertices as the strengths of the entanglement constraints on a length scale 27r/fc 
on the tracer dynamics. The time dependence described by the normalized propagators 
captures the disentanglement processes. The net friction is obtained from the summation of 
friction amplitudes on all length scales weighted by the characteristic rates required for decay 



of the corresponding constraints^. Let us first recall the previous results where the matrix 
was assumed to be dynamically frozen in the PMC model^'^^'^^. This corresponds to setting 
Sk{t) = Sk in eqn. (15). The relaxation of the entanglement friction then only proceeds 
via the probe collective dynamics calculated in the RR model. In the asymptotic limit of 
large degree of polymerization of the tracer the intrinsic (A^-independent) characteristic 
length scales like 6, ^p, and a become widely separated from Rg oc y/N — > oo. Thus, eqn. 
(15) further simplifies as the wavevector dependence of the matrix structure factor, Sk, and 
the direct correlation function, c^, can be neglected. Also, the asymptotic result of eqns. 
(34) and (35) can be used, and the PRISM Ornstein-Zernicke equation, Sk — cok + Qmhk 
together with eqn. (13) leads to the replacement^^ S^QmCo ^ —1. The suppression of the 
center-of-mass diffusion coefficient, Eq in ©qn. (16), then is of the form^^'^^'^^: 



The entanglement friction is proportional to the square of the number of binary segmental 
contacts of a pair of interpenetrating polymer chains weighted by a factor proportional 
to the probability of contact, i.e. the intermolecular pair correlation function at excluded 
volume distance. In PMC theory the fluctuating forces describing the entanglement friction 
therefore arise from excluded volume interactions of strongly interacting polymer chains. 
Their strength is given by the number of two macromolecules pair-contacts, approximately 
N'^/Rg, and the probability of close contact of two monomers on different macromolecules, 
Qd- The terminal relaxation time and viscosity follow as <x rj oc Rg/ D oc g'^N^ . 

It is well known from neutron or light scattering experiments that the collective density 
fluctuations of a polymeric melt rapidly decay into equilibrium^^"^^. At least within the 
accuracy of the mentioned experiments, Sk{t) has decayed to zero at times much shorter 
than To ■ Naively, this contradicts the frozen matrix assumption in eqn. (15), which requires 
that Sk{t) is nonzero on time scales where the tracer dynamics determines the decay of the 
entanglement friction^. However, obviously it is not necessary that the full amplitude of the 
collective density fluctuations is frozen in. Strongly molecular weight dependent contribu- 
tions to the tracer dynamics are already obtained if some small, but finite, A^-independent 
amplitude of the matrix dynamics only relaxes at times of the order of the probe dynamics 
as given by eqn. (33). Indeed, this is the experimental situation for the stress relaxation 




(43) 



function G{t), where Gn <^ QmkBT. Little is known experimentally and theoretically about 
the existence, magnitude and time dependence of an entanglement plateau in Skit) in the 
melt^^^^^. In theta solutions, however, the plateau in Sk{t) has been observed and there the 
general expectation has been verified that its amplitude is proportional to the ratio of shear 
to compressibility modulus^^. Note that for nearly incompressible melts the sensitivity of 
neutron or light scattering is not sufficient for direct observation of such small amplitudes. 

The question of the entanglement amplitude in the slow matrix dynamics motivated 
the re-derivation of the PMC equations projecting onto the collective stress variables in 
section 2. A. Experimentally, the entanglement plateau in the stress modulus is well studied^. 
In this section it will be shown that straightforward assumptions lead to identical PMC 
expressions starting from either the choice of matrix density or matrix stress fluctuations as 
a slow collective variable, and that the previous frozen matrix assumption nicely connects 
to established concepts. 

Let us first comment on our treatment of binary polymer-solvent solutions which we 
crudely model as an effective one-component polymer ffuid. Adding a small molecule liq- 
uid to a polymeric liquid has two kinds of effects. First, the equilibrium structure of the 
polymeric subsystem is changed. The polymer density fluctuations grow and the osmotic 
compressibility strongly increases^'^'^^'"^^'^^. Locally, the contact probability of segments on 
different macromolecules, i.e. gd, goes down^^'^^. Also the intramolecular structure can be 
changed depending on the type of solvent^'^^'^^. We will use a simple Gaussian chain model 
and an integral equation approach, PRISM^^"^^, in order to capture these effects as well 
as possible and study their consequences on the dynamics of entangled polymers^^^^^. A 
second group of effects of adding a solvent directly influences the dynamics of the polymer 
component. The local mobility, i.e. is strongly affected^. The two component liquid 
mixture exhibits further transport processes, interdiffusion^^ or possibly gel-like modes^*^'^^, 
which may couple to the polymer density ffuctuations. Hydrodynamic interactions arise 
from the instantaneous but long ranged solvent motions^'^'^^. The local mobility affects 
our overall time scale and is taken from comparison with experiments. The remaining two 
effects described above are neglected. The standard argument is that the interdiffusion is 
fast following the solvent motion and that hydrodynamic interactions are screened^. Our 



approach therefore focuses on the effects of the hquid equihbrium structure on the dynamics 
but does not modify the basic dynamical equations with varying solvent content and quality. 

The long lived matrix constraints can be written as a wavcvcctor dependent amplitude 
and a normalized time and (in general) wavevector dependent function: 

Skit) = SkftMt) , where « r,) = 1 . (44) 

The normalized correlation function, ^k(t), describes the final decay of the matrix con- 
straints which possess the frozen-in amplitudes Skfk ■ Faster decay processes, the so-called 
microscopic and glassy relaxation processes, are assumed to be completed. They result in a 
decay of Sk{t) from its initial value Sk{t — 0) — Sk down to the arrested amplitude, S^f^ 
where f§ < 1. The importance of an amplitude of the final relaxation step smaller than 
unity has most clearly been recognized in the mode coupling theory of the glass transition^^. 
There equations for the frozen in amplitudes are derived. In the present context we will 
describe a physically plausible model and quote neutron scattering measurements in order 
to obtain a simple model for f^. 

The amplitude of the matrix constraints contains two spatially varying or wavevector 
dependent factors which describe what might be called the "porosity" of the entanglement 
constraints. First, the equilibrium matrix structure, Sk, describes a nonuniform compress- 
ibility of the surrounding polymeric liquid. We assume a simple Ornstein-Zernicke form 
which is appropriate for concentrated or (semi-) dilute theta solutions and for wavelengths 
large compared to monomeric length scales^'^'^^'^^. 

where is the density screening length or "mesh-size", and Sq = QruksTuT, where ht is 
the isothermal (osmotic) compressibility of the polymer subsystem. The mesh size can be 
measured by neutron or other scattering techniques^. Alternatively, liquid state theories 
predict its magnitude and density dependence in agreement with blob scaling and field 
theoretic considerations^^'^^'^^. In PRISM theory the form of eqn. (45) and the additional 
result Cfc ~ Co has been derived^^'^^. Whereas a melt is described by a density screening 
length of the magnitude of monomeric length scales, a less concentrated polymeric liquid 
exhibits a much larger ^p. grows with decreasing concentration until, in the dilute solution 



limit, it agrees with the radius of gyration of a single polymer chain^. In good solvents the 
correct form of the density susceptibility nontrivially deviates from the Ornstein-Zernicke 
form^ of eqn. (45). As we require Sk in order to perform wavevector integrals only, we 
neglect these quantitative corrections. In this extended blob-scaling picture of de Gennes 
the difference of a good and a theta solvent is only captured in the different dependence of 
monomeric length scales on the actual polymeric density. Whereas in a theta solvent the 
statistical segment length scales are density independent, in a good solvent the excluded 
volume interaction leads to an effective monomeric size dependent on density^'^^, a ~ Q^^^- 
In eqn. (44), only an amplitude /jf < 1 of the total density fluctuations is arrested up to 
the disentanglement time . We follow Semenov^^ and Genz^^ in estimating the magnitude 
of its spatially homogeneous part, /,f , as the ratio of the shear modulus to the bulk modulus, 
or in our notation = So/Ne- Note that in typical melts where Sq ~ 0.25 and N^, = 50 — 
300 one finds J'q is of the order of 10"^ or smaller. We shall not repeat the arguments of 
Semenov and Genz, but rather derive the same result for Zl(t) looking at the shear stress 
variables below. Clearly, the very local density correlations should not be arrested by the 
entanglement effect^'^^. We follow the general conviction, as it is captured for example in 
the reptation/ tube model or in the PMC theory, that there exists an entanglement length b, 
6^ = NeCr"^, which determines the spatial resolution of the matrix entanglement constraints. 
Adapting the results of Ronca^^, de Gennes^^ and des Cloiseaux^^ about the plateau in the 
collective single chain dynamic structure factor, to the plateau in the collective matrix 
structure factor, we arrive at the following simple model: 

/f = {Som e-(^^/6)^ . (46) 

The entanglement length in the mentioned theoretical considerations and in the comparison 
with neutron scattering is defined by fk=Q/b = We have checked explicitly that the 

simplified Gaussian shape assumption in eqn. (46) does not affect the results of our theory 
appreciably. 

The sketch in figure ^ of the matrix and tracer correlations and characteristic length 
scales summarizes the physical factors which enter the PMC memory functions. Clearly, in 
order to obtain the N ^ oo asymptotic results the tracer size, Rg, has to greatly exceed 
the matrix correlation lengths, C,p of the density fiuctuations and b of the elastic mesh. It 



will be a central part of our results and their discussion that the intramolecular correlations 
of the tracer coarse grain over a much larger matrix volume in the conformational friction 
function, M(t), than in the ccntcr-of-mass memory function, S(t). 

The form of the time dependence of the matrix constraints in Sk{t) can be deduced fol- 
lowing the theoretical studies mentioned above^^'^^'^^'^^. The matrix provides entanglement 
constraints up to a time when the final disentanglement step of the matrix polymers takes 
place. The reptation/ tube picture for the collective single chain correlator, and the PMC 
results for the generalized Rouse mode correlators^, indicate that this disentanglement pro- 
cess is characterized by a uniform, wavevector and Rouse mode independent relaxation time 
To . The natural assumption for which also agrees with the theoretical findings of 

Semenov^^ and Genz^^, therefore is^: 



Using the eqns. (44-47) in eqn. (15), and the original expression for the M-memory 
function^'^^'^^ (containing Sk{t)), one obtains the final expressions of the PMC theory for 
the tracer diffusion coefficient, the disentanglement time To , and the shear viscosity of 
linear chain polymers^'^^ . Since the disentanglement time is determined from the M- 
memory function, eqn. (30), but also enters in its relaxation, via eqns. (27) and (47), a 
self-consistency aspect emerges. Before proceeding to its discussion and consequences for 
transport properties of the PMC theory, it is instructive to derive the same final equations 
based on the assumption that the collective stress variables are slow, i.e. starting from eqn. 



Eqn. (27) was obtained from assuming that the collective matrix stress variables are the 
slow variables hindering and slowing down the tracer dynamics. It is well known that the 
final stress relaxation is characterized by the disentanglement time scale r^, so that we can 
write^"^ with eqn. (47): 



Neglecting small quantitative differences, the glassy elastic modulus can be approximated 



^^(t) = $(i) = e-(*/^^) . 



(47) 



(27). 



(48) 



by: 



or 



G, 



'k ~ QmkBTiSk/So) 



QmkBT 



(49) 



Eqn. (49) expresses the idea that the glassy, elastic stresses are of the order of ksT per 
segment and their correlation length is the density screening length; note, that this expression 
is only needed for inverse wavevectors large compared to the monomeric sizes. It is well 
established experimentally that the entanglement plateau in the shear modulus is a factor 
l/Nf. smaller than the glassy modulus^. We again follow the generally held idea that the 
entanglement length b determines the spatial correlations of the elastic constraints provided 
by the entanglements. A simple model for the wavevector dependent plateau in the shear 
modulus therefore is^^~^^: 



It does not appear to be a coincidence that the original PMC equations considering the 
collective density fluctuations and, the "new" PMC equations resulting from eqns. (27) and 
(47) — (50) which consider the collective stress variables as slow, lead to identical expressions 
for the memory functions, E(i) and M{t), and consequently to the same theoretical results. 
Rather, we surmise that apart from small quantitative differences identical expressions for 
the fluctuating force memory functions are obtained if the same physical assumptions are 
made. Namely, (1) the fluctuating intermolecular forces are projected onto the product 
of a tracer and a collective dynamical variable. (2) The resulting four point correlation 
function is factorized as in eqn. (24), stating that the friction forces relax via the collective 
matrix dynamics and the collective tracer density fluctuations. (3) The collective tracer 
dynamic structure factor is evaluated within the RR model. (4) There is some long-lived, 
finite amplitude of the collective matrix constraints which only relaxes during the final 
disentanglement process. Its spatial correlation length is b, the entanglement length. (5) 
The disentanglement time follows from self-consistently requiring the internal relaxation 
time to agree with the collective flow time. 

Two important new ingredients result from the more detailed treatment of the matrix 
constraints as compared to the previous frozen matrix considerations. First, the constraint 
release mechanism, i.e. the consideration of the time dependence of the matrix constraints, 
relaxes the unphysical assumption that some collective degrees of freedom do not relax into 
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(50) 



equilibrium. Instead, the known slow dynamics of the disentanglement process enters the 
theory in a self-consistent way. Moreover, the variation of the tracer dynamics with matrix 
molecular weight can now be studied. Second, the treatment of the constraint porosity is 
improved. Whereas some previous numerical studies^°'^^ included the spatial variation of the 
compressibility via 7^ 0, the frozen-in elastic mesh, characterized by the entanglement 
length 6, enters as a new length into the PMC fluctuating force memory functions. An 
immediate consequence of 6 7^ is the prediction of much larger finite size corrections since 
the inequality b ^ holds in general, and finite size corrections within the PMC theory 
only vanish for Rg large compared to all other length scales^'^^. 

Note, that the entanglement degree of polymerization, A^e, can be estimated from the 
theory, i.e. from the crossover calculation of the RR model, eqn. (35). However, it can 
also be obtained experimentally from a rheological measurement; = QmkBT /N^.. In the 
following we will use the prediction of the present version PMC theory that entanglement 
corrections enter as a function of N and in a scaled fashion N/Ne, section 4.B discusses 
the reasons for this. Thus, we assume to be known experimentally and express results 
in terms of N/N^. 

4. Asymptotic and General Finite Size Results 
A. Asymptotic Predictions 

The results of the PMC model in the asymptotic, N ^ 00, limit have been worked out 
previously^'^^'^^'^^. We first recall these results for the transport coefficients. 

Changes of the polymer transport coefficients depending on tracer molecular weight 
result from the long-ranged contributions {k oc Rg^) of the entanglement friction. The 
matrix constraint amplitudes and their friction contributions in the S- and M-memory 
functions differ due to the different intramolecular weighing factors, k'^LUk and k'^uil in E(t) 
and M{t) respectively. For Gaussian chain polymers a simple but rather accurate expression 
for the intramolecular structure factor is^: 



Figure ^ shows the intramolecular factors of the entanglement amplitudes in the S- and 
M-memory function. Global entanglement constraints, on the size of the tracer, dominate 
in the conformational memory function, M{t), only. The diffusive tracer dynamics, eqn. 
(33), enhances the small wavevector contributions to the Markovian, transport coefficients. 
Figure ^ shows that for conformational dynamics and stress relaxation long wavelength 
[kRg ^ 1) contributions extending across the size of the tracer molecule contribute most. 
In the center-of-mass friction local modes dominate the amplitude associated with 
entanglement constraints, although kRg ^ 1 contributions dominate the long time (Markov) 
friction coefficient renormalization. In PMC theory, this difference in weighting of force 
correlations is the origin of the different, property-specific molecular weight scaling of the 
corrections to the Rouse, single polymer dynamics which is so characteristic of entangled 
polymeric systems. 

The diffusion coefficient follows from the uniform drag memory function, S(t) eqns. (15) 
and (16) plus the new analysis of section 3.B: 

D = [1 + Xn^]-\ (52) 

where the Rouse diffusion constant, = [ksT / C^qN), is the result for a single macro- 
molecule immersed in a continuum system characterized by an instantaneous friction 
coefficient'^'^^, Co- A reptation like asymptotic scaling is obtained as the corrections to the 
Rouse result are of quadratic order in the interaction parameter^^'^^'^^ (gdVN). The strength 
parameter of the PMC corrections relative to the Rouse result, i.e. Xc, is proportional to a 
well defined, and in principle independently measurable, quantity: the mean squared, aver- 
aged intermolecular force exerted by the matrix polymers on the probe center-of-mass per 
unit density: 

32 

Ad = - ~ {\F\')/g„^ • (53) 

It is important to note that in general or a, which is defined in eqn. (53), are expected 
to be density and material dependent parameters. This result differs from the reptation 
prediction where Ad = 3 is a universal number independent of solvent character, chemical 
structure, and density of the polymeric system'^. 



The inverse strength parameter a will play the most important role in quantifying the 
finite size corrections in the PMC approach, and can be more explicitly written as^'^^: 

- = ^{QMS,f ~ {m/Qm . (54) 

Note that in previous PMC work, e.g. ref. 33 and earher pubhcations^'^'^^"^^'^^, the hteral 
frozen matrix approximation lead to an expression for a of the the form: l/a oc g^N^. 
As will be discussed in the following paper no significant changes in previous results for the 
concentration dependence of polymeric transport coefficients arise from this difference based 
on PRISM input for simple Gaussian chain models. The actual magnitudes, however, of the 
PMC dynamical parameters now come out in almost quantitative agreement (perhaps for- 
tuitously) with experimental measurements. The parameter a also determines the prefactor 
of the asymptotic final disentanglement time, and by assumption, of the shear viscosity: 

^ = ^^(1 + A.(^)^), (55) 



where the Rouse result, -rf- — QmCo'^^N /2>Q, again follows for a single Gaussian polymer in 
a continuum with friction coefficient Co but without hydrodynamic interactions^. Here is 
connected to a via: 

A, = - , (56) 
a 

which again can be compared to the reptation prediction^ of A^ = 15/4. These equations 
predict a universal ratio of the diffusion constant and shear viscosity in the asymptotic limit: 

^ - Gn^ - gA for TV ^ oo . (57) 
tig oAd Id 

Remarkably, the value of the constant X^i/GXd = 3/16 is very close to the reptation result 
of 5/24. Note, however, that the prefactors Xd and A^ are in general material and ther- 
modynamic state dependent, and that their close connection via the parameter a follows 
from two key points of the PMC theory. First, the treatment of the translation-rotation 
coupling discussed in context with eqns. (11) and (26), and second, the assumption that the 
entanglement friction in the uniform and the conformational memory functions, S(t) and 
M{t), decays via the identical tracer dynamical process, i.e. the RR collective intramolecular 
density fiuctuations. Of course, the neglect of hydrodynamic interactions and the simplified 



treatment of solutions possibly could also affect eqn. (57). The exact numerical factors in 
eqns. (53) and (56) depend on the exact form of the intramolecular form factor Uk used in 
the vertices. The values denoted above result when, for analytical convenience, the simple 
approximation eqn. (51) with Rg — is used. 

The asymptotic result for the conformational relaxation time, r^, , is directly connected 
with eqn. (55): Tu/tq — N^/{aNe). This result may also be taken as the (experimental) 
definition of the (inverse) entanglement strength parameter a. Note, that the reptation/ tube 
result. To — 3r^(iV/iVe), differs from this PMC result in the prediction of an universal value 
of the asymptotic prefactor, a — tt^/S for reptation, whereas a is given by eqn. (54) in PMC 
theory. Dielectric relaxation, which measures the end-to-end-vector correlation function, 
{P{t) • P(0)), does not exactly determine , as it weights global modes much more heavily 
than viscoelastic measurements. From the spectrum of modes obtained in PMC theory after 
neglecting the tracer chain ends as discussed in section 2. A and refs. 6, 31, and 34, one can 
easily calculate an averaged dielectric relaxation time, — dt{F{t) ■ P(0)) / (P(0) • P(0)) . 
In a dielectric spectrum, it describes the low frequency behavior, e"{u!) oc cut^ for cut^ <C 1. 
We find 

t'/t" = ^[1 + Kn] , where = 32/q; . (58) 

The difference in the asymptotic hmit, t^/to — > 32/12, arises from the contributions of 
the low lying modes to the end-to-end-vector fluctuations which are described by S'(i) in 
eqns. (1) and (26). Naturally, the approximation to neglect the chain end effects when 
diagonalizing eqn. (1) will lead to the largest systematic errors for the low lying modes and 
therefore in quantities like r^. 

B. Parameters and Magnitudes of Finite Size Corrections 

The asymptotic predictions for the transport coefficients of a tracer derived above assumed 
frozen-in constraints of the surrounding polymer matrix. A simple additivity approximation 
was used in order to describe the crossover from unentangled Rouse to entangled PMC 
dynamics. In the rest of this manuscript and the following paper, we will discuss corrections 
to these results arising from the spatial correlations and the time dependence of the matrix 



constraints or entanglements. The aim is not to describe more accurately the low molecular 
weight regime, < 5A^e, but to study the effects which prevail up to the highest molecular 
weights which are experimentally accessible. As will be discussed, finite corrections arc 
predicted, which vanish for A?" ^ oo, but, depending on the situation studied, are not yet 
neghgible for N/N,, 10^. It is for such corrections at rather high N/N^. > 5, that we 
hope our simple models for the matrix constraints and the Rouse to entangled dynamics 
crossover to be appropriate. As well defined equilibrium structural information is required, 
improvements can easily be considered if better (experimental) information is available. The 
deviations of the transport coefiicients from the Rouse values depend on molecular weight via 
the ratios n — N/Ne and p — P/Pe only, where N (Ne) is the tracer and P (Pe) the matrix 
polymer degree of polymerization (entanglement degree of polymerization). The physical 
origin of this scaling is the existence of the entanglement length, b. Entanglements form an 
elastic mesh of spacing b and slow down the polymer dynamics only if the intramolecular 
correlations of the tracer extend over distances larger than b. The entanglement friction of a 
polymer of size Rg therefore depends on the ratio ~ N/Ng. This result, however, rests 

upon our corrected identification of the frozen-in amplitudes of the entanglement constraints, 
i.e ~ So/N^ in cqn. (46) or ~ l/N^ in cqn. (50). 

The asymptotic results, eqns. (52) and (55), are determined solely by the rescaled 
molecular weight, n, and the intermolecular excluded volume strength factor, 1/a. From 
the comparison with the reptation/ tube predictions for the diffusion constant and viscosity 
it appears that a value oi a ~ 3 should be considered reasonable^ at least under melt 
conditions. In the following papcr'^^ we will give further arguments using PRISM results 
and comparisons with experiments for this choice. The finite size corrections depend on a 
second, dimensionless material dependent parameter, S, where: 

6 = Ub . (59) 

S is the ratio of the compressibihty length scale (physical mesh), ^p, to the elastic mesh size, 
b. The appearance of this ratio again follows from the spatial origin of entanglements within 
the PMC theory. The constraints arise due to equilibrium, eqn. (45), and dynamically 
frozen-in, eqn. (50), intermolecular correlations. Again, we generally expect 5 is a material 
and density dependent parameter. As will be discussed in the following paper^^, from the 



recent comprehensive neutron scattering and rheological study and theoretical PRISM 
and computer simulation results, values of and h are known for a variety of polymers 
and varying density conditions. Based on this data, we estimate 5 ~ 0.05 for melts and 
somewhat larger values for solutions^^, 5 < 0.3. 

PMC theory differs from the reptation picture in that it generally predicts a material 
dependence in the asymptotic transport coefficients even when expressed in reduced vari- 
ables. Moreover, the appearance of the structural correlations, parameterized by S, has no 
analog in the tube description, which sets effectively 6 = 0. This statement does not concern 
the scaling considerations for the asymptotic transport coefficients nor plateau moduli^'^^'^^, 
but addresses the finite size corrections, i.e. the speeding up of the dynamics, due to the 
correlations of the entanglement constraints. 

Introducing the reduced variable, q = kRg, the two memory functions of the PMC 
approach, E(t) from eqn. (15), and M{t) from eqn. (27), can be written using the two 
non-universal parameters a and 5: 



= /^(g) S{qS/V^) f{q/V^) e-t^^l e'^*/--) , (60) 

To A*^ ^0 



M{t) = r dq I'^iq) S{qd/V^) f{q/V^) e^^^^' e'^*/"^) , (61) 

where the reduced compressibility, S = Sk/So, follows from eqn. (45), and the elastic 
frozen-in amplitude, / normalized to /(O) = 1, from eqn. (50). The first time dependent 
factor in eqns. (60) and (61) arises from the tracer, eqn. (33), the second from the matrix 
dynamics, eqn. (47). The amplitudes of the entanglement forces are proportional to the 
spatial weights: 

= 7^ • ^''^ 
^"(«' - Mf. • 

where the rescaled intramolecular structure factor, cug = a;(fc=g/^g)/A^, appears. An im- 
portant quantity specifying the relevance of the constraint release mechanism compared to 



the tracer dynamics is the ratio, denoted as ©, of the tracer RR time scale to the matrix 
disentanglement time: 

e(n,p,a,S)= (64) 

where — / dqq^Cj'^S{q5 / ^/n) / J dqq^u^ describes the finite size corrections in the 

RR Markovian calculation, eqn. (34); 7(3;) — > 1 for x — > 0. Eqn. (64) indicates one source 
of variations if tracers of chemistry different from the matrix polymers are studied. The 
monomeric friction coefficients and the appropriate entanglement molecular weights will 
differ in general. However, also the vertex and consequently a, eqn. (54), will differ, as 
the total intermolecular site-site correlation function, hk eqn. (13), is changed by different 
tracer-matrix coupling depending on the chemical interactions^^. Therefore, we restrict our 
considerations to chemically identical tracer and matrix polymers in the following ( except 
for section 6); i.e. = 1 in eqn. (64). 

The factor /?(p, a, S) in eqn. (64) describes the suppression of the disentanglement time 
relative to its asymptotic value: 

r./ro = , (65) 

where /3 — >• 1 for p — >• 00 follows from eqns. (30) and (61). (3 is found as the solution of the 
self-consistency equation requiring that the matrix disentanglement time agrees with the 
longest internal relaxation time of one of the matrix polymers. 

Pip, a, S) = r dq /^(g) S{qS/Vp) /(g/VP) J^^^fl ^2 fr, " ^^^^ 

If the matrix disentanglement time, , is found from eqns. (64) to (66), the same equa- 
tions determine the tracer internal relaxation time with the obvious replacements, especially 
Q{p,p, a, 5) — > Q{n,p, a, S) in eqn. (66) and p ^ n elsewhere. 

The shear viscosity including the finite size corrections and the crossover to the unen- 
tangled Rouse result then follows as: 

77 = 77^[l + A,p2/3(p,a,5)] . (67) 

Note that our crossover model qualitatively, but presumably not quantitatively, correctly 
describes the viscosity at low molecular weights. The longest relaxation time, Tu , however. 



does not crossover to the Rouse time, — tq{P/t^Y, with decreasing P, but goes to zero 
below a certain P. It is at such low molecular weights, P ^ Pg, where the self-consistency 
equation, eqns. (64) and (66), loses its validity. 

From the so-determined matrix disentanglement time, , the diffusion coefficient of a 
tracer can be calculated using eqns. (16) and (60): 

D^^[l + Xnn K{n, 5, e{n,p, a, S)) , (68) 

where the reduction of the asymptotic entanglement friction constant follows from the E 
memory function, eqn. (60): 



K{n, 6,Q)= dq I^{q) S{q6/V^) f{q/V^) — ''^ ' ^ . (69) 

Jo S{n,p,a,d) + q^/cuq 

Again, the crossover to the Rouse result is described qualitatively, but possibly not quanti- 
tatively, correctly 

In order to compare with dielectric loss measurements of the end-to-end-vector 
fluctuations^^ we must include the finite size corrections into eqn. (58) for the relaxation 
time, T^, of the end-to-end-vector correlation function of a tracer in a polymer matrix. 
Two qualitatively different finite size corrections appear in as the conformational mem- 
ory function, M{t), but also the homogeneous friction function, S'(t), affect the low ly- 
ing, global modes of an entangled polymer. It is shown in refs. 6, 31, and 34 that the 
generahzed Rouse mode correlators of PMC at zero frequency have the following form: 
C'^(uj = 0)/Cp{t = 0) = To[{^f{l + E(j) + Mo]. One notices the mode independent relax- 
ation rate, TqMo which dominates the viscoelastic dynamics, and the Rouse-like rates, which 
show the familiar mode dependence, and are negligible except for the lowest modes. 
From this, summing over the odd modes^, p — 1,3..., of the end-to-end-vector correlation 
function at zero frequency, one obtains the dielectric relaxation time, r^: 

n 

T^ = To — [l + (24/?(n,p,a) + 8 K'{n,p,a,S) ) - ] , (70) 
12 a 

where K' is the normalized finite size correction arising from the time integral over the 
uniform friction function, E'(i): 



K' differs from K{n, p, a, 6), eqn. (69), in the contributions from tlie intramolecular structure 
only: 

The quantity /3 in eqn. (70) follows from eqns. (64) and (66) with the replacements denoted 
there, and the numerical prefactors in eqn. (70) assure the correct limiting behavior, eqn. 
(58), as /3 — >• 1 and i^' — > 1 for p — > oo. Again, we suggest eqn. (70) for intermediate and 
large molecular weights of the tracer, n, and of the matrix polymers, p. When comparing 
to experimental data it has to be kept in mind that our present results apply for chemically 
identical tracer and matrix polymers only. System specific variations will arise from differ- 
ences in the monomeric friction coefficients, in the molecular weights of entanglement ( see 
eqn. (64)) and from different local packings, i.e. ga, and effective potentials, Cq, in eqn. (54). 

It is one of the central results of the PMC theory that the entanglement effects on 
the center-of-mass motion and the conformational dynamics are closely connected via the 
two memory functions^'^^''^^ S(t) and M{t). They differ in the intramolecular weighting 
factors, If/N"^ eqn. (60), and I^^ eqn. (61), only. Note from figure ^ that , eqn. 
(63), is peaked around q = kRg ^ 1, whereas J^, eqn. (62), monotonically increases with 
increasing q. Friction contributions on the global, macromolecular size scale determine the 
internal relaxations. This is one of the central findings of PMC^, which also agrees with the 
reptation/ tube idea that a motion of the whole chain is necessary in order to relax internal 
fiuctuations^''^. This behavior of J*^ strongly differs from the (tightly connected) finding 
that contributions to the friction of the center of mass motion, i.e. eqn. (62), arise from 
more local correlations^^. See figure ^ where the entanglement amplitudes and the friction 
contributions, i.e. the amplitudes weighted by the corresponding relaxation rates (in the 
asymptotic, frozen matrix limit), are shown. 

This difference causes three major consequences of PMC theory. First, the asymptotic 
slowing down of the internal modes is larger by a factor A^^ than the extra friction exerted on 
the center-of-mass^. Second, the porosity or spatial correlations of the matrix constraints, 
Sk < Sq and < f^, influence the center-of-mass motion more strongly. The uniform 
drag friction function, S(t), is much more sensitive to the actual values of the flnite, ma- 
terial dependent length scales, and b. The conformational friction function, M{t), on 



the other hand is relatively independent of these non-universal parameters. Third, the en- 
tanglement friction can decay via the constraint release mechanism, i.e. the decay of the 
matrix constraints following eqn. (47), or via the tracer collective structural dynamics of 
the RR model, which exhibits a diffusive pole, eqn. (33). Comparing the two decay rates 
for different wavevectors, the tracer dynamics will dominate the decay at large wavevec- 
tors, q = kRg ^ 1, and the matrix djTiamics will dominate in the limit q = kRg 0. 
Therefore, the constraint release mechanism will lead to much larger finite size corrections 
in the conformational friction function, M(t), than in S(t), as a much larger portion of the 
entanglement constraints arises for q = kRg < 1 in M than in E. The internal dynamics and 
the viscosities, therefore, will be dominated by the constraint release mechanism in general, 
whereas the constraint porosity will mainly influence the center-of-mass motion. In the 
following sections we discuss these differences and other aspects of the finite size effects in 
more detail by numerically solving eqns. (60) — (69) for the physically relevant parameter 
ranges, a ~ 3 and 6 =0 to 0.3. 

5. Predictions for Transport Properties and Relaxation Times 
A. Effects of Constraint Release on the Viscosity 

Considering a melt of polymers with degree of polymerization N = nN^, the non-linear 
equation for the finite suppression, j3{n,a,6) eqn. (66), of the disentanglement time can 
be easily solved numerically. (3 depends on the molecular weight via the dimensionless ratio 
n = N/N(. only. The non-universal parameters a, where 1/a ~ {\F\^) / Qm in eqn. (54), and 
5, eqn. (59), lead to in general density and chemical structure dependent results. Results 
for (3 are shown in figure |^ spanning a wide range of possible variations in a and especially 
5. Since in melts b ~ 35 — 90A depending on chemical structure''^, it is apparent from 
figure ^ that the density screening length has to be increased to unphysically large values, 
5 = ^p/b ~ 0.5, in order to observe changes in j3 in the range n > 5, where our equations 
are reliable. From the discussion in section 4.B of which spatial correlations dominate the 
respective friction functions, it can easily be understood that the internal relaxation time 
To is only affected by the constraint release mechanism but not by the constraint porosity. 



The disentanglement time, Tn , and as a consequence the viscosity, rj, follow from the 
conformational friction function, M{t), which is most sensitive to constraints for q = kRg ^ 
1. These rather long range, intermolecular correlations decay more effectively via the matrix 
dynamics than via the diffusive RR tracer dynamics. The local spatial correlations of the 
entanglements, either the equilibrium correlation length, ^p, or the localization length of the 
frozen-in amplitude, b, influence the conformational dynamics only very weakly since the 
local contributions to M{t) are virtually negligible. 

Note that for physically reasonable choices''^'^^ of S, i.e. for 6 < 0.3, almost no effects on 
P or To can be seen in figure ^ for n > 5. For n < 5 smaller molecular weights our crossover 
model cannot be considered very reliable. Thus, we conclude that for experimentally relevant 
parameter ranges the disentanglement time, r^, , or viscosity, r], is suppressed below its 
asymptotic value by the constraint release mechanism. To a good approximation the self- 
consistency equation, eqns. (64) and (66), can therefore be simplified by neglecting the 
constraint porosity, i.e. setting Sq = 1 and = 1 in eqn. (66). The correction factor (3 
then depends on the ratio n'^/a only, i.e. l3{n,a) = (3{n/a'^). Its asymptotic behavior can 
be found easily: 

A 2 

f5{n/a'') 1 - \{—Y'^ for n oo . (73) 

The solid line in figure ^ indicates the solution, /5(n/a^), to the self-consistency equa- 
tions neglecting constraint porosity. The asymptotic behavior, eqn. (73), and the in- 
dependence on 5 for large n can be observed when comparing the solutions, /3(n, a,(5), 
of eqns. (64) and (66) including constraint porosity. A handy interpolation approxima- 
tion to (3{n/a^) in the relevant molecular weight range is given by the following formula: 
f3{n/a^) = exp-[x + 0.22x^ + OMx^ + 0.02x^] where x = 4/3 (^)^/^ in figure | it lies on 
top of the exact curve, (3{n/a^). In the following, i.e. after figure ^ in this paper and in 
the following paper, we will simplify the numerical work and neglect the constraint porosity 
corrections in the conformational dynamics. The above simple but accurate expression for 
(3{n/a'^) will be used. Errors are thereby made in regions only where the simple matrix and 
crossover models we employ are not reliable anyway. 

Figure |^ shows the predicted viscosities corresponding to these disentanglement times 
where the thick lines are given by the simplified expression neglecting the spatial constraint 



correlations: 

r/ = r/^[l + \n^p{n/a^)] . (74) 

Again, in figure § it is seen that only unrealistically large values of 5 lead to appreciable 
deviations of r] from the constraint release result, eqn. (74). Also, even though eqn. (74) 
does not lead to a rigorous power law behavior except in the asymptotic limit, r] ~ A^^ for 
A^ — s> oo, an effective power law over two orders of magnitude with an exponent of ~ 3.4 can 
be accurately fitted to the results. Note that the extremely slow approach to the reptation- 
like asymptote, rj ~ A^^(l — c/n^^^)^ where c = (4/3) -^a, leads to considerable deviations 
even for n ^ 10^, molecular weights which are difficult to achieve experimentally^^. In figure 
^, for values n = 10^, the correction factor, j3, is still of the order 0.7 or less. Experimentally, 
of course, the very slow drift of (3 in that range may not be discernible in the data and an 
apparent power law with an exponent approaching 3 may be concluded for n > 10^ — 10'^. 
Note that the apparent power law rj ~ A^^'^ holds in a molecular weight range, n > 10, 
where our crossover model and the considerations about the matrix structural functions 
may be oversimplified, but, are qualitatively correct, and contain all the physics expected 
to arise in full numerical PMC calculations. As our prediction is not a rigorous power law 
in the intermediate n range, it is difficult to estimate what effective exponents could be 
reported when fitting power laws to eqn. (74). The results will also depend on the window 
in molecular weight of the fits. Exponents up to 3.5 seem easily achievable. 

Two effects of the inverse strength parameter a on the viscosity can be seen arising from 
the constraint release mechanism. First, asymptotically the viscosity becomes proportional 
to 1/a. Obviously, the stronger the intermolecular excluded volume forces are, the higher 
the viscosities will be in the entangled regime. Second, increasing a leads to more curvature 
and higher effective exponents in t] as the PMC corrections to the Rouse result effectively set 
in at larger n only but then increase more rapidly with n. Finally, since our predicted finite 
A^-corrections to the viscosity are essentially independent of 6 = ip/h, we predict a sort of 
universality of the noriasymptotic behavior nearly independent of both chemical structure 
and polymer concentration when r] is expressed in terms of N/Ng. This prediction appears 
to be in excellent agreement with experiments^'^^. For the values of n = 10 — 1000 relevant 
to experiments, we find that is typically in the range < 1 — 0.1. Thus we conclude that 



there is no wide separation between the characteristic single chain conformational relaxation 
time and the time scale for entaglement force decay. Hence, a non-Markovian situation is 
suggested under experimentally accessible conditions. 

Let us comment on another, non-Markovian consequence of the constraint release mecha- 
nism, which will affect the PMC results for the frequency dependence of the shear modulus at 
rather low frequencies, cuTd < 1. Neglecting the matrix dynamics, i.e. the constraint release 
relaxation, eqn. (47), and setting $(t) = 1 in eqn. (61), the conformational friction function 
exhibits a long time taiP'^^'^'', M{t > r^^) \/N {T^'^'/tf/^, where r«« = ToN^y/n. This 
slow decay leads (for finite A^) to an anomalous, low frequency behavior in the storage part 
of M{u), M'{uj) ~ v^(a;T««)=^/^ which is more weakly dependent on frequency than the 
M' ~ cu'^ expected from a Markovian separation of time scales. In numerical solutions to 
the full PMC equations^^, neglecting the constraint release, it has been observed that this 
anomalous decay of the conformational friction function leads to a final disentanglement 
peak in G"{uj), which is quantitatively broadened for ujto < 1 compared to the experimen- 
tal data^°~^^. Obviously, the exponential matrix constraint relaxation, eqn. (47), yields 
M{t) oc e"*^*/"^^) (see eqn. (61)) which will cutoff the long time tail in M at the disentangle- 
ment time, To . Therefore, we expect the quantitative discrepancy^^ between PMC results 
and experiments for shear moduli at low frequencies, ujTo < 1, will be resolved in future 
numerical solutions of PMC theory including the constraint release mechanism. Note, that 
it also can be expected that no changes will be found for the high frequency power law 
tails of the disentanglement process, loTo ^ 1, where the analytic PMC prediction^^'^^'^^ of 
exponents around 0.2 — 0.25 and the numerical PMC results'^^'^^ agree nicely with shear 
measurements^°^^^ but differ strongly with the pure reptation model prediction of a;~^/^. 

B. Constraint Porosity Effect on Diffusion Constants for P ^ N 

Upon immersing a tracer polymer with degree of polymerization N — nN^ into a polymeric 
matrix of degree of polymerization P = pN^, the constraint release mechanism can be 
neglected if p ^ n. Then the matrix polymers are immobile relative to the tracer. In eqn. 
(64), the ratio of the corresponding time scales, ©, vanishes as is evident from the factor 
(n/p)^. The entanglement friction in the memory functions then can relax only via the tracer 



collective structure factor (computed with the RR model). The resulting tracer diffusion 
coefficients will in general still show finite n corrections as the matrix constraints are not 
homogeneously correlated, i.e. h/Rg and ^p/Rg are nonzero. The tracer center-of-mass 
experiences the full entanglement constraints only if it has to distort the frozen-in elastic 
mesh appreciably {Rg ^ h) , and to compress the equilibrium structure uniformly {Rg ^ ^p) . 
As the friction contributions to the conformational dynamics arise from q = kRg ^ 1, 
the internal conformational degrees of freedom are rather insensitive to the spatial finite 
size correlations. Therefore, in the limit P ^ N ^ Ng, the PMC approach predicts the 
reptation like scaling, r^/ro = N"^/ (aNe), for the internal relaxation time of a tracer polymer. 
The center-of-mass motion, however, because the more local intermolecular correlations 
contribute heavily, eqn. (62), does not feel the full entanglement constraints as long as the 
finite length scales, ^p, b are not negligible relative to Rg. The diffusion constants therefore 
are increased relative to the N ^ oo asymptotic behavior. Moreover, even at fixed scaled 
degree of polymerization N/Ng corresponding to Rg/b"^, we find that the density screening 
length, although small, leads to a non-universal reduction of the friction. Note that the 
screening length, (at least in semidilute solutions), and the entanglement length, b, are 
much larger than monomeric sizes, i.e. cr or a persistence length for flexible polymers. 
Therefore, parameterization of the dynamical consequences of these length scales with a few 
chemistry dependent parameters, as done in our models eqns. (45) and (46), seems justified. 

Figure |^ shows the tracer diffusion constants, D^, for two values of the mean square 
force per density parameter, eqn. (54), a = 3 and a = 1/2, corresponding to Xd = 32/9 
and A/) = 64/3 respectively. Values of the screening to entanglement length ratio, 6 = ^p/b, 
are chosen in a wide range. In the first case of A/? ~ 3.6, which we suggest to be a melt 
like case, a very rapid crossover from the Rouse, D ~ 1/A^, to the asymptotic PMC result, 
D ~ I/^XdN"^) from eqn. (52), is seen for realistic melt values of 5 < 0.05, Importantly, 
in the range n > 5 almost no deviation from the asymptote can be observed. Increasing S 
beyond its physically expected melt range, somewhat larger exponents in an effective power 
law, reaching up to D ~ A^~^'^ may be observed in the intermediate n-regime. In the second 
case, Xo ~ 21, which we consider relevant for solutions (see the following paper), where also 
6 < 0.3 is expected, an intermediate behavior appears in the range 1 < n < 50. There 



the diffusion constant lies appreciably above its asymptote. For example, for Xo = 64/3 
and 6 = 0.3 {6 = 0.05) at ri = 5 the diffusivity is increased by a factor DX^n/D^ 3.7 
(1.8), respectively. In this n-range, a free fit to the numerical results would also lead to 
A^-scaling exponents exceeding the classical reptation result appreciably. In figure |^ a power 
D ~ N""^'^ is drawn for comparison. The crossover to the true asymptote happens at much 
larger molecular weights in the solution case than in the melt case. 

Note that these deviations from the reptation like scaling have no analog within the 
extended reptation/ tube ideas. Constraint release is irrelevant because of the arrested 
matrix dynamics, P oo, and the contour length fluctuation mechanism of Doi^^'^^ does 
not affect the A^-scaling of the tracer diffusion constants. It appears that a microscopic 
approach incorporating the liquid structure properly is required to identify these finite size 
corrections arising from the spatial correlations of the matrix entanglement constraints. 
Chemistry and density or composition dependence is expected for these effects. Also note 
that in contrast with the reptation prediction, the asymptotic prefactor of the PMC, Xd, is 
non-universal and is expected to be polymer-density dependent. 

C. Variation of Tracer DifFusivities with Matrix Molecular Weight 

The variation of the tracer diffusion constant upon changing the matrix molecular weight can 
also be described within our general approach. Because of the neglect of the matrix dynamics 
in the RR calculation the theory is restricted to the range p = P/Pe > 1. First, the matrix 
disentanglement time has to be found from the self-consistency equation, eqns. (64) and 
(66). Then the suppression of the entanglement effects on the center-of-mass motion of the 
tracer can be calculated from eqns. (68) and (69). Figure |^ shows the case with parameters 
argued to be relevant for polymer melts. With increasing P a rather rapid crossover of the 
tracer diffusion constant from its Rouse value, D^, to the strongly entangled value, 
discussed in section 3.D, is observed. The self diffusion coefficients, for N = P, are also 
indicated and generally fall in the transition region between the two asymptotes, and 
D^^. In figure H, the parameters a = 3, 6 = 0.05 and = 3, lead to a little overshooting 
of the asymptotic relative to the rept at ion-like asymptotes. 

The dependence on matrix molecular weight in the crossover region is not a rigorous 



power law but, if approximated by one, corresponds to exponents of the order of 2. This 
exponent is smaller than the limiting ideal behavior of D ~ P^^ expected from Grassley's 
constraint release formulation^^ and D ~ expected from Klein's analysis^^'^^. However, 
if different material parameters are chosen, in figure ^ we use a = 3, 6 = 0.3 and = 18 
corresponding to dense solutions, then larger effective exponents are found. Clearly, the 
dependence on the matrix molecular weight must be weaker than the one following from 
D^"^ ~ l/Tr,{P), as this result would hold only if the tracer dynamics was totally frozen in 
and solely the constraint release mechanism, or matrix dynamics eqn. (47), would relax the 
entanglement constraints in eqn. (60). In this solution-like case a stronger matrix 

molecular weight dependence in the transition range is observed compared to the melt case. 
D ~ p-2.8 jg (^Yamn for comparison in figure |^. Again, the self diffusion coefficients lie in 
the crossover region. For the case of figure the constraint porosity is also important and 
leads to deviations of D^^ from the reptation like predictions for all matrix molecular weights 
as discussed in section 5.B. Note that the values of the diffusivities in the solution case lie 
below the corresponding (same n = N/N^.) melt diffusion constants because of the larger 
prefactor, A^, of the large asymptote. 

In figure 5 (and slightly also in figure H), for the sake of discussion the theoretically 
predicted correlation of A/j with a, eqn. (53), was violated. In this case a determines the 
matrix disentanglement time r^, and the viscosity via the eqns. (61) and (64) — (67) but 
not the strength factor. Ad, in the S memory function, eqn. (60), describing the tracer. 
If the matrix viscosities corresponding to figures ^ and ^ were plotted, they would agree, 
but the tracer internal relaxation time would differ from the matrix one. A experimentally 
relevant situation where such a breaking of the relation in eqn. (53) can appear is in the 
case of a tracer which differs chemically from the matrix polymers. Then the intramolecular 
correlation function, hk in eqn. (13), which determines the vertices or constraint amplitudes, 
is expected to be different. That is, intermolecular tracer-matrix packing is not the same as 
the pure matrix-matrix correlations. Specific tracer-matrix-polymer chemical interactions 
caused by the addition of selective solvent could also lead to such a more complicated 
situation. 



D. Variation of Tracer Dielectric Times with Matrix Molecular Weight 

The dynamics of a tracer polymer in a polymeric matrix is slowed down by the entanglement 
constraints. In the limit of rapidly moving matrix polymers, the Rouse model with instanta- 
neous friction coefficient, Co, would apply for the tracer. Dielectric spectroscopy^°'^^ measures 
the slowing down of the relaxation time of the end-to-end-vector correlation function, 
of eqn. (70). Asymptotically the reptation like scaling, r"^ ~ A^^, eqn. (58), is predicted by 
PMC theory. 

As the global modes contribute to appreciably, in general r*^ differs from the internal, 
conformational disentanglement time . Whenever it is justified to neglect the contribu- 
tions from the uniform drag friction, i.e. ck' = %K'nja <^ cp = 24:(3n/a in eqn. (70), 
then the relation = 2r£, holds. In this case, the conformational friction determines the 
dielectric relaxation time, , and therefore results independent of the spatial constraint 
correlations, constraint porosity, but strongly dependent on the constraint release mecha- 
nism are predicted by PMC theory. Whenever the uniform friction arising from S'(t), eqn. 
(26), dominates the dielectric relaxation time, strong constraint porosity effects resulting 
from the finite lengths scales, b/Rg > and C,p/ Rg > 0, are predicted in the PMC approach. 
Note that the neglect of chain end effects, discussed in section 2. A, strongly affects the 
global mode contributions in eqn. (70). Without exact numerical diagonalizations of the 
full equations of motion, eqn. (1), no estimate of the resulting errors in can be given. As 
figure |Ty shows this uncertainty mainly affects the results for in the limit of large tracers 
in entangled matrices, n ^ p ^ 1. Melt-like parameters, a = 3 and 6 = 0.05, are used for 
figure 

Our model, eqn. (70), rests upon the assumption of strictly identical matrix and tracer 
polymers, which may differ in degree of polymerization only. Due to the approximations in 
the RR model the results apply for entangled polymeric matrices, p = P/Pe > 1 only. Also, 
the curves are physically relevant for r^j/r-^ ^ ^ 1 only, where the connection between 
the conformational friction function and the final time scale, eqn. (30), holds. Note that 
r^, from eqn. (65) does not, as would be required physically, crossover to the Rouse time 
in the limit N ^ P but vanishes asymptotically because of cp ~ p^/n^^"^ for n ^ p. This 
limitation is not a fundamental one, but rather a consequence of our simplified modeling of 



the Rouse to entangled crossover. 



In figure [10| one notices that apart from small quantitative corrections the conformational 
friction determines the dielectric relaxation time, , for tracers in strongly entangled ma- 
trices, i.e. for n < p. For various matrix molecular weights, the dielectric, , and internal, 

, relaxation time follow the asymptotic, r/r-^ ~ n, behavior for n <^ p. For n > p, the 
uniform friction contributions to dominate as they increase asymptotically, ck' ~ p^/^n^/^ 
for n ^ p ^ 1. For long tracers in shorter, but still entangled matrix polymers, PMC theory 
therefore predicts an intermediate power law behavior, ~ A^i^/s for n ^ p ^ 1, resulting 
from the contributions of the global modes to the dielectric spectrum. The matrix dielectric 
relaxation times agree with the r*^ of a tracer of identical degree of polymerization, n = p, 
and show a r^, ~ A^^'^ behavior for large n, as is also observed in the viscosity for this choice 
of strength parameter, a = 3. 

The constraint release mechanism and the constraint porosity thus strongly affect the 
molecular weight scaling of the dielectric relaxation time r^. Note, however, that the high 
frequency wing of the disentanglement process, as described by PMC theory, is connected 
to the tracer shape fluctuations. The physical origin of the fractal frequency behavior, e.g. 
(P"(u;) ■P(0)) ~ io^^/^, eqn. (42), or (P"(c<j) ■P(0)) ~ cu"^/^^ for higher frequencies, therefore 
is predicted to be very different from the one determining the r"^ scaling with N. In the 
limit P ^ N, where reptation like scaling in the time scale, r ~ N^, is found since the 
matrix constraints cannot decay appreciably, the shape fluctuations of the tracer still speed 
up the early stages of the disentanglement process relative to pure reptation and lead to a 
much reduced high frequency slope for P"{uj) at intermediate frequencies, c^r^ ^ 1; see the 
discussion in section 3. A. 



6. Extension of PMC Theory to Polymer Transport through Random, Gel Like Media 

The diffusion of polymer tracers through amorphous media is of considerable interest in 
both physical and biological science On the one hand, the original reptation theory 
was formulated for the situation of a polymer in a gel^. On the other hand, gel electrophore- 
sis constitutes a powerful but still poorly understood technique to characterize polymers, 
mostly biological ones''^. In the present section we cannot address the multitude of different 



realizations of amorphous materials or of finite driving fields in electrophoresis. However, the 
unique aspect of PMC theory is that it connects the dynamics of polymers to the underlying 
structure of the surrounding matrix. This allows one to study a number of interesting tracer- 
gel cases within the same formalism. Moreover, the structural information, i.e. the spatially 
resolved compressibility and the mesh of elastic constraints, that uniquely determine the 
transport properties of a polymer tracer (if its intramolecular correlations are specified), 
raises the possibility of predicting the tracer dynamics from purely static experimental data. 

In PMC theory an amorphous solid or gel is characterized by the time independence 
of the matrix constraints the probe feels. That is, in eqn. (44) or eqn. (48), $(i) = 1. 
Crosslinked polymeric gels or silica gels exhibit arrested density structures which do not 
relax into homogeneous, liquid-like equilibrium. The constraint release mechanism therefore 
is not operative. 

In order to describe gels which are fractal and exhibit pores on all (relevant) length scales 
a simple fractal ansatz for the density correlations is used: 

^ (1 + A;2^2^Pf/2) ' C^^) 

where is a characteristic mesh or pore size. Operationally the gel fractal dimension is 
defined by a scattering experiment in the intermediate, self-similar spatial range, k^g S> 1, 
where Sk ~ /c"^^ is seen. Exponents in the range Dp 5/3 — 2 apply to polymeric gels 
prepared by crosslinking (good) polymer solutions^. 

Note that eqn. (75) describes a very simple model of a gel and in reality may not 
be adequate for experimental systems. Effects hke quenched disorder of crosslinks^" are 
not described by eqn. (75), but can be put into Sk if the appropriate experimental or 
theoretical information is available. It should be mentioned that quenched crosslinking 
disorder introduces heterogenieties^° on an even larger length scale than which results in 
enhanced values of Sk at small k. For example, the following form has been shown to fit 
experimental data on some gels^° 

1 + k^eg 

where ^xc > Cg- the spirit of our approach, which does not differentiate between quenched 
and annealed disorder, such a structure in Sk would lead to further nonuniversal finite size 



effects in the polymer tracer dynamics extending to even higher tracer molecular weights as 
found in the present study. 

The density screening length, ^g, can be considered to be the size of the largest pores in 
the gel. It separates the homogeneous density structure from the self-similar intermediate 
range. Of course, the breakdown of eqn. (75) at microscopic length scales is neglected for 
the dynamics of entangled polymers. 

Two distinct classes of gels can be defined following the discussion of the constraint 
porosity in section 3.2. First, in a structurally rigid or "hard" gel, the mass density structure 
can be thought of to pose strong constraints on the motion of a probe. The full amplitude of 
the gel density fiuctuations therefore is expected to contribute to the entanglement friction 
functions. This formally corresponds to: 



in eqn. (44). Silica gels presumably belong to this class of hard gels^^. A second class includes 
"soft" gels prepared from crosslinked polymer solutions. Even though the entanglement 
constraints then are permanent due to the crosslinking reaction, it cannot be expected that 
the strength of the so-formed elastic mesh significantly exceeds the strength of the non- 
crosslinked, time-dependent precursor. The model of section 3.2 should therefore apply and 
the constraint amplitudes again contain a small factor: 



For lack of detailed knowledge about the spatial correlations of the elastic mesh in crosslinked 
gels, we continue using a Gaussian ansatz for the wavevector dependence of f^. Different 
/c-dependent functions will affect the comparison with experiments only in so far as shifts 
in the poorly a priori known (fit) parameters occur. 

In the soft gel models, the ratio ^g/b for crosslinked polymer gels is expected to be of the 
order of the values found for non-crosslinked polymer solutions^ {S — ^p/h ~ 0.3). Larger 
values are expected in the so-called strangulation regime, where the high crosslinking density 
forces the entanglement length to become smaller than the one caused by the temporary 
crosshnks^^ . 




(76) 



(77) 



The question of the intramolecular structure of a polymer in a fractal medium is very 
complex and still not well understood^^. Ideal, Gaussian-hke intramolecular correlations 
cannot be expected in general. The mass-size scaling exponent, u, where Rg ~ A^'^, may 
even depend on degree of polymerization, i/ — i^{N), and on other system-specific features^®. 
We will neglect these difficulties, and assume a fractal model with fixed mass-size exponent, 
v, for the tracer polymer as well: 

where c — {1 + 2i/)(l + z/)/ci/(2r(l + l/2i/))^^. From the obvious generalization of the 
E memory function, eqn. (15), one obtains the following results for the tracer diffusion 
coefficient in a fractal gel in d-space dimensions: 

D'' = ^[ 1 + N'-''''Kg{Rg/^g, Ub) , (79) 



where 



/ dqq'^-'u;'^S{qUR ,)f'{Qb/R,) I dqq<^-^u;'^S{qCj Rg) 
(/ dqq'^-^ujlf 



Kg{Rg/ig, ig/h) = ^ ^ , ,.^d-lr,2^2 ^ " (^O) 



The dynamics is predicted to depend strongly on the ratio Rg/^g. This dependence of the 
tracer diffusion coefficient on the ratio of the tracer size relative to a mesoscopic length scale 
which characterizes the matrix mesh, is closely connected to the constraint porosity effects 
discussed in section 5.2. The latter have been shown to be important in polymer solutions 
or if special tracer-matrix chemical interactions exist. 

The strength parameter Xg in eqn. (79) can, in principle, be determined from the tracer- 
gel structural correlations; i.e. the pair contact value and the direct correlation function, 
Cfe. For example, the PRISM description of binary polymeric liquids can be evaluated in 
the limit of vanishing tracer concentration^^'^^. However, because of the possibly varying 
chemical interactions of the tracer with the gel pore walls, no general statement about the 
magnitude of is possible. Moreover, the increase in the amplitude f^, eqn. (77), which 
likely results from the crosslinking, also increases Xg, but cannot be described theoretically 
at present. Note that in the true polymer solution case, i.e. chemically identical tracer and 
matrix polymers, eqn. (77) applies, v — 1/ Dp — 1/2, and Xg is connected to Ad, eqn. (53), 
via Xg — Xo/Ne- For more rigid gels, however, larger strength parameters, Xg, are expected. 



Moreover, all A's are proportional to the tracer-gel interaction strengths and therefore are 
expected to increase if the tracer polymers adsorb to the gel pore walls, or if other specific 
tracer-gel or tracer-solvent interactions exist. Much larger values than Xo/Ne therefore are 
possible for A^. 

From eqn. (79) the asymptotic scaling of the tracer diffusion coefficients with molecular 
weight can be inferred, as the finite size correction factor, K, eqn. (80), approaches unity for 
large tracers, Rg ^ C,g- Asymptotically, D ~ ]\f-5+2diy jg fQ^^d independent of the gel fractal 
dimension^^''^^. In 3-dimensions for ideal coil tracers the classic D ~ A^~^ law is recovered. 
Note that eqn. (79), in the same way as our melt and solution descriptions, includes a 
crossover model describing the tracer diffusion constant for all tracer molecular weights; 
for the unhindered tracer the Rouse model is again chosen to apply. The PMC results 
including finite size effects therefore present a unified description of the tracer dynamics in 
gels arising from the competition of free and entangled polymer motion, where the second 
contribution depends on the length scale ratio Rg/^g. In this respect the PMC theory 
differs from phenomenological approaches like the entropic barrier model of Muthukumar 
and Baumgartner^^'^^^*^^, where special geometric considerations for Rg ^ are invoked, 
and the asymptotic limits {Rg ^ ^g or Rg -C ig) and smooth regime crossovers are not 
included in the description. 

Figure |TT] shows tracer diffusion coefficients resulting from eqns. (79) and (80) for fractal 
dimensions corresponding to self avoiding random walk polymers, v = l/Dp = 3/5 in the 
Flory approximation^'^'^. Gel pore sizes, ^g, are denoted by stating the degree of polymer- 
ization of a tracer whose radius of gyration agrees with C,g, i.e. C,g = Rg{N). Results for 
two different values for are shown in order to explore the variation with this unknown 
system-specific parameter. 

In figure ^ one notices that deviations from the Rouse, unconstrained diffusion start 
for Rg/ig below 1 but do not lead into the asymptote, Kg{x ^ oo,y) ^ 1, up to very high 
degrees of polymerization of the tracer. Intermediate effective power law behavior arises 
with varying exponents, D ~ N'"^'^ to D ~ A^~^'^, depending on the parameters varied in 



figure 0. The steepest molecular weight dependences fall into a range where the tracer Rg 
exceeds the gel pore size ^g appreciably. This finding appears to disagree with the entropic 



trapping ideas where maximal non-asymptotic results are centered around Rg ~ ^g. 
In the case of PMC theory the finite size effects lead to anomalous exponents for much larger 
tracer sizes. 

The results of eqns. (79) and (80) very weakly depend on the fractal dimension of the 
gel. Assuming for simplicity the identical constraint strengths, i.e. in eqn. (79), the soft 
gel, eqn. (77), and the hard gel, eqn. (76), lead to similar intermediate effective power laws. 
Somewhat higher exponents result for soft gels. Increases in the strength parameter, Ag, 
lead to steeper crossover curves. Note that the existence of the entanglement length in the 
flexible mesh systems or soft gels leads to extended Rouse like behavior as the tracer has to 
uniformly distort the elastic mesh, and therefore has to be larger, in order to feel the full 
entanglement constraints. 

The degree of polymerization dependence of the diffusion coefficients, for intermediate 
tracer sizes, cannot be rigorously described by power laws. If, however, power laws are fitted 
to the numerical results, effective exponents clearly exceeding the reptation value of 2 are 
obtained. Note that this holds although the true asymptotic molecular weight dependence, 
D ~ A^~^/^ in figure is even weaker than the reptation-like scaling. An upper bound 
to the effective exponents can be deduced from the small tracer size limit in eqn. (80). 
For Rg/ig one easily finds Kg ~ {Rg/^gY^^ ~ jyavD^^ which in the limit of a large 
asymptotic prefactor, A^, and or large pore sizes leads to the upper bound, D ~ ]\j-^+'ii^{d-Dp) 



where d is the spatial dimension. In the case of figure 11 this estimate is ~ A^' 



-3.4 



7. Discussion 



The predictions of PMC theory for the transport properties of entangled polymer melts, 
solutions, and gels are determined by the intramolecular and intermolecular equilibrium 
structural correlations. A number of approximations are necessary in order to simplify the 
A^ generalized Langevin equations of eqn. (1). 

1. The fiuctuating intermolecular forces are approximated by their statistical overlap with 
the collective tracer structure factor and a collective matrix correlator in accord with known 
mode-coupling ideas. It is expected that for any choice of slow matrix variables qualitatively 
identical PMC expressions result, as shown for the examples of matrix density and stress 



fluctuations. 

2. The slow dynamics of the intermolecular forces is connected to the diffusive dynamics 
of the tracer and the (generaUy nondiffusive) disentanglement process of the surrounding 
matrix polymers. 

3. The (projected) tracer dynamics entering the PMC friction functions is taken from a 
short time and/ or small molecular weight calculation, the RR model. 

4. The matrix or constraint release dynamics is evaluated from a self-consistency argument, 
requiring that the single polymer and matrix disentanglement time agree. 

5. Simple models for the strength of the entanglement constraints exerted by the polymer 
matrix are formulated. They lead to qualitatively identical results when considering various 
slow matrix variables, as shown for the examples of density and stress fluctuations. 

6. In order to calculate bulk transport properties or response functions for polymeric liquids, 
it is assumed that the entanglement effects arise from incoherently added, single chain 
contributions. For example, the bulk viscosities and dielectric susceptibility are derived 
from a single chain calculation. 

7. Solutions are treated in the same way as melts, neglecting special dynamical effects, like 
"gel-modcs"^'''®^, but concentrating on the equilibrium structural changes. 

Assumptions 6. and 7. arc familiar from the reptation/ tube theory, and are generally 
motivated by the peculiarity of the entanglement effects to macromolecules (point 6.), and 
the similarity of these phenomena in melts and solutions (point 7.). 

The robustness of the PMC results with regard to the identiflcation of the slow matrix 
variables constraining the tracer polymer (points 1. and 5.) rather reassuringly verifies that 
the structural correlations of the matrix constraints are characterized by two length scales, 
the entanglement length, b, and the density screening length, ^p. Whereas the first length 
scale is postulated and identified as the tube diameter in the phenomenological reptation 
approach, the effects of the second length scale on the dynamics are neglected there. In 
PMC theory, the matrix constraints are not fully developed if the ratio 6 = ^p/b is not 
small. Smaller effective entanglement friction and finite size corrections then arise and 
vanish for large molecular weights only. The concept of "lack of full topological correlation of 
entanglement constraints" due to Muthukumar and Baumgartner^^ bears some resemblance 



to our ideas, but there are significant differences. 

The idea to equate the single chain and collective disentanglement times (point 4.) ap- 
pears rather natural and is also suggested by point 6. In a crude sense it is related to 
the constraint release mechanisms of Grassley^^ and Klein^^'^^ within the phenomenological 
tube model framework although there are strong differences. For example, these tube-based 
approaches are not self-consistent, and significantly affect only tracer diffusion and not 
conformational relaxation or viscosity in the N = P case. 

Employing the RR collective dynamical structure factor (point 3.) in the PMC friction 
functions is suggested by the break down of the Rouse model for degrees of polymerization 
above A^e- Actually, A^e is defined by calculating where the Rouse model as a zeroth order ap- 
proximation is overwhelmed by the first order correction, the RR model^^. For the molecular 
weights and times required in the PMC friction functions, i.e. for N and t 3> r'^lNe), 

the RR correlator is physically more reasonable than the Rouse correlator. Nevertheless, 
this aspect of the theory seems the most "uncontrolled" . The use of a fully self-consistent 
correlator, i.e. entering the PMC correlator into the friction functions of PMC theory, is 
not appropriate, if naively implemented, since it results in a severe overestimation of the 
entanglement friction and an arrest ( "macromolecular scale localization transition") at a 
finite degree of polymerization^^'^^. 

The determination of the decay of the intermolecular forces from the collective dynamic 
structure factor of the tracer bears close connection to one of the central assumptions of 
reptation: the coherent center-of-mass motion of the chain in the tube determines the con- 
formational relaxation^''^. In agreement with this picture, PMC theory finds, when neglecting 
tracer shape fluctuations, identical results for the conformational dynamics as reptation. Not 
only the A^-scaling of the internal relaxation time and viscosity agree, but also the shear 
stress spectrum shows the high frequency asymptote, G"{ujTu ^ 1) ~ o;"^/^, familiar from 
the tube survival function of reptation^'^. In better agreement with measurements of shear 
and dielectric spectra, PMC results^^'^^'^^ including tracer shape fluctuations find more shal- 
low slopes, G"{ujTo 3> 1) ~ u}~^ , where x ^ 0.2 — 0.25. In agreement with recent dielectric 
measurements of Adachi and coworkers^°'^^ , these non-reptation like power law behaviors 
are predicted to persist even, if, for strongly entangled polymer matrices, the reptation like 



scaling of the internal relaxation time, ~ N^, is observed. 

Whereas the tracer shape fluctuations only influence the initial decay of the disentan- 
glement process, there arise two finite size effects in PMC theory which affect the transport 
coefficients of entangled polymers. Both are closely connected as they arise from the con- 
sistent model for the matrix entanglements discussed in section 3.B. First, a decrease of 
the entanglement friction results from the time-dependent decay of the matrix constraints, 
termed constraint release mechanism. Clearly, this idea bears similarity to Grassley's^^ or 
Klein's^^'^"^ ideas for tracer diffusion (but not viscosity), although the actual implementa- 
tions strongly differ. PMC theory describes the decrease of the effective friction coefficient 
whereas in refs. 15 — 17 two independent relaxation rates (or diffusion constants) are added 
based on the assumed statistical independence of the reptation and constraint release trans- 
port processes. Second, the spatial correlations of the matrix constraints, termed constraint 
porosity, enhance the tracer motion as the full constraints are effective only if the inter- 
molecular forces are summed over finite spatial regions. The constraint correlations are of 
two distinct origins: spatial compressibility, characterized by the density screening length 
or mesh size, C,p, and elastic mesh or entanglement length, b, correlations. 

One of the central findings of the PMC approach is the result that entanglement con- 
straint contributions to the conformational dynamics arise from globally, across the tracer, 
correlated intermolecular forces, whereas the center-of-mass friction results from more local 
intermolecular forces®. The corresponding intramolecular factors determining the constraint 
amplitudes are either peaked around kRg = 1 or increase monotonically, see figure ^. This 
finding appears connected to the reptation/ tube idea that the mobility of the polymer in 
the tube can be motivated by the picture of "pulling a wet rope through a tube"^. The mo- 
tion of the whole "rope" is slowed down by short ranged fricition. Whereas, conformational 
relaxation requires the decay of the tube correlations due to the motion of the chain ends^'^. 
Also, the aspect of an underlying diffusive collective tracer dynamics, i.e. the collective RR 
dynamical structure factor, agrees with this picture. 

A first consequence of the difference in the spatial correlations of the conformational, 
M{t), and the uniform, S(t), memory functions is their very different Markovian value®, 
Mq ~ and Sq ~ A^. Because of this feature of strongly different effects of the entan- 



glements on the internal and the center-of-mass friction, the prediction of conformational 
and stress relaxation arrest (plateau) but continued (but slowed down anomalous) segmental 
diffusion follows. 

A second consequence, worked out in this paper, is the very different sensitivity to finite 
size effects predicted by PMC theory for internal and center-of-mass dynamics. Again, 
because of the rather long ranged spatial correlations of the entanglement friction in the 
conformational dynamics, constraint porosity is irrelevant there. However, the constraint 
release mechanism strongly speeds up the conformational dynamics as it overwhelms the 
slow tracer RR dynamics on long length scales or for small wavevectors. The center-of-mass 
motion, on the other hand, is accelerated due to the decrease of the entanglement friction on 
local, finite length scales. The constraint porosity, therefore, enhances the diffusion constants 
of tracer polymers even in matrices of effectively immobile matrix polymers. 

Two nonuniversal, microscopically defined parameters, the inverse entanglement strength 
parameter, a, and the length scale ratio, 6 = ^p/b, are predicted to control the finite size 
corrections to the ratios of the transport coefficients relative to their Rouse values if molecular 
weights are normalized by the corresponding entanglement molecular weight. The first 
parameter exists in the reptation/ tube approach, where an universal value, a ~ 3-4 is 
calculated. The second parameter has no analog there. 

PMC theory predicts an extremely slow approach of the disentanglement time, r^, , 
and hence the viscosity, to their asymptotic, ~ r^, ~ A^^, behavior. Higher exponents 
result from fits of effective power laws to the numerical PMC results for molecular weights 
n = M/Mf, < 10^. No clear separation of the final disentanglement time from the time scale 
for entanglement force decay is expected in experimentally relevant parameter ranges. 

Diffusion constants for polymer tracers in matrices of immobile polymers are unaffected 
by the constraint release process but differ from the asymptotic behavior, D ~ A^~^, because 
of the spatial variation of the entanglement constraint amplitudes. Figure |^ exemplifies 
that for the experimentally relevant intermediate molecular weights effective exponents, 
D ~ A^~^, with X significantly exceeding the classical value of 2 will often occur. 

The variation of the tracer diffusion constants upon varying the molecular weight of 
the matrix polymers results from a combination of both types of finite size effects. Strong 



variations are obtained due to the constraint release mechanism for systems where the matrix 
polymer molecular weight is smaller or slightly larger than the tracer molecular weight. 

Constraint release and porosity also compete in the finite size corrections of the end- 
to-end-vector relaxation time, , measured in dielectric spectroscopy^^'^^ . In the limit 
of immobile matrix polymers, the disentanglement and dielectric times are proportional 
and follow reptation-like scaling, r^, ~ ~ N^. In the limit of large tracer polymers in 
entangled matrices, a new behavior of is predicted, ~ N^^^^. It results from the global 
mode contributions in the dielectric spectrum and requires M ^ P ^ Mg, where M, P and 
Me are the tracer, matrix polymer and entanglement molecular weight, respectively. 

The dynamics of tracer polymers in amorphous solids, especially in gels, does not require 
quahtatively new physical effects to be appended to the PMC description once the physically 
appropriate generalizations of the matrix and tracer structure are made. Obviously, in an 
amorphous solid the entanglement constraints cannot fully relax, and hence the constraint 
release mechanism is absent. The constraint porosity, however, leads to strongly molecular 
weight dependent tracer diffusivities, especially, as the tracer-gel interactions may often lead 
to stronger intermolecular forces than in simple homopolymer systems. For non-Gaussian 
intramolecular correlations of the tracer polymer, Rg ~ N'^, PMC theory does not predict 
reptation like results, but rather D ~ _/V-5+2d;/ asymptotically in d spatial dimensions. Very 
much stronger intermediate molecular weight dependences are predicted, however, arising 
from the constraint porosity corrections as long as the ratio of gel pore size, ^g, to tracer 
size, Rg, is not negligible. In the limit of large tracer-gel interactions an upper bound, 
D ~ ]\j'~5+2v{d-DF) ^ may be approached, where Dp is the fractal dimension of the gel. These 
stronger effective molecular weight dependences do not indicate in any obvious manner 
the existence of qualitatively different mechanisms of polymer transport, but emerge as 
natural generahzations of the constraint porosity effects predicted by PMC theory also for 
polymer melts and solutions. Indeed, this is a very important point applicable to all our 
new results. Traditionally, a change in A^-scaling exponents is interpreted in terms of a new 
dominant transport mechanism at the level of individual polymer trajectories. In contrast, 
PMC theory focuses on the entanglement friction, and the diversity of possible scahng laws 
emerge as a consequence of the influence of structure and multiple competing length scales 



and relaxation channels on fluctuating force time correlations. Real space physical motions 
do not need to be guessed as in phenomenological theories, nor can we unambiguously infer 
them. 

In the following paper these theoretical predictions will be tested by quantitative compar- 
isons with experimental data. A necessary preliminary is the estimation of the equilibrium 
parameters specifying the dynamics which will be accomplished using integral equation the- 
ories, computer simulation results, and experimental scattering and rheological data. 
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Figure Captions 

• Figure 1: Schematic representation of the (static) entanglement constraints in the 
uniform-drag friction function, T,(t) eqn. (15), of PMC theory. The tracer, char- 
acterized by the intramolecular correlations co and length scales a and Rg, interacts 
via the short-ranged (excluded-volume diameter d) pseudo-potential, c (the direct 
correlation function) with the matrix polymers. The entanglement constraints of the 
matrix separate into local (density screening length ^p) compressibility and mesoscopic 
(entanglement length b) elastic correlations. 

• Figure 2: Schematic figure showing the restriction of the segment dynamics for 
times shorter than the disentanglement time. A discrete bead-spring polymer model 
is used, where eqn. (31), with the unknown Cq, set to = for simplicity, allows 
the determination of all bead displacements from the motions, u(t) and v(t), of two 
arbitrary segments. 

• Figure 3: Schematic representation of the (static) entanglement constraints in the 
friction functions of the PMC theory. The tracer interacts via the short ranged pseudo- 
potential c (the direct correlation function) with the matrix polymers. The entan- 
glement constraints of the matrix separate into local (density screening length ^p) 
compressibility (structure factor S) and mesoscopic (entanglement length b) elastic 
correlations (amplitude /). Note that only a small number of the matrix polymers 
are shown as the blobs of size which fill space. The tracer exhibits different in- 
tramolecular correlations, /, for the uniform drag, ~ cu, and for the conformational 
friction, ~ uj'^/k'^. Note that the conformational friction coarse grains the matrix 
correlations over much larger spatial regions. 

• Figure 4: Intramolecular factors determining the entanglement constraint amplitudes 
in the conformational, (long dashes) eqn. (63), and in the center-of-mass, (solid 
line, scaled by 1/10) eqn. (62), friction functions plotted versus reduced wavevector, 
q — kRg. The corresponding asymptotic friction amplitudes, P^ = 2I^uj/q^ (chain 
curve) and P^ — 2I^uj/q^ (short dashes), follow from the diffusive relaxation rates, 
2Cjq/q^, of the RR model in eqn. (33). 



Figure 5: Finite size reduction factors, P{n,a,6) = ToaNe/roN^, of eqn. (65), for 
different entanglement strengths, {\F'^\) / ~ and length scale ratios, 6 = ip/h. 
The thin lines correspond to the parameters a as denoted and to increasing parameter 

5 from left to right. The solid line shows the solution, (3{n/a'^), neglecting constraint 
porosity, and the dotted line the asymptote eqn. (73). 

Figure 6: Ratios of shear viscosity to Rouse viscosity corresponding to the disentan- 
glement times To following from the results of figure |[ Bold lines denote the results 
neglecting constraint porosity for different a and the thin lines for fixed parameter a 
shift with increasing parameter 5 from left to right. Lines of the same style belong to 
one value of parameter a. A power law, rj/r]^ ~ A^^ '', is shown for comparison. The 
N ^ oo asymptote \r,{N/NeY is drawn as a dotted line for a = 2 or = 6. 

Figure 7: Tracer diffusion constants in the limit of immobile matrix polymers, P 3> 
A^. Two sets of curves are shown with two different values of the asymptotic prefactor, 
\di and corresponding to upper or lower horizontal scale. The curves show increasing 
steepness in the intermediate region with increasing 5 = ip/h; 5 = 0.01, 0.05, 0.1, 0.2, 
0.3 and 0.5 respectively. Three power laws with arbitrary prefactors are shown for 
comparison. 

Figure 8: Tracer diffusion coefficients in a polymer melt as a function of matrix molec- 
ular weight, P/Pe, for different tracer degrees of polymerization, N/Ne] logio{N/Ne) = 
0, 0.4, 0.8, 1.2, 1.6 and 2 from top to bottom. The parameters employed, a = 3, 

6 = 0.05 and Xd = 3, are argued to describe a polymer melt. The self diffusion 
constants, V, are denoted by circles and long dashes, and can be compared to the 
asymptotic, D^'^ for P ^ oo, values shown with horizontal short dashes. A power 
D ~ P^^ is shown for comparison. 

Figure 9: Tracer diffusion coefficients in a polymer solution as a function of ma- 
trix molecular weight, P/Pe, for different tracer degrees of polymerization, N/Ne] 
logio{N / Ne) = 0, 0.4, 0.8, 1.2, 1.6 and 2 from top to bottom. The parameters em- 
ployed, a = 3, 6 = 0.3 and Ad = 18, are argued to describe a polymer solution or a 
melt exhibiting a special tracer, matrix polymer interaction. The self diffusion con- 



stants, (circles and long dashes), the asymptotic values, D^^ for P — > oo (short 
dashes), and a power D ~ P~^-^, are shown. 

Figure 10: Tracer dielectric relaxation time, r^, normahzed by the Rouse time as a 
function of reduced tracer molecular weight, n = N/N^, for different matrix molecular 
weights, p = P/Pe, p = S, 10, 30, 100, 300, and 1000 from bottom to top. Melt like 
parameters, a — 3 and S = 0.05, are chosen. Chemically identical tracer and matrix 
polymers are considered. The asymptotes, ~ A^^ for n <S p (long dashes) and 
~ N^^/^ for n » p (short dashes), are included. The matrix end-to-end-vector 
relaxation times, shown with circles and a dashed-dotted line, correspond to n = p. 
A dotted line shows a power law, ~ A^^'^, fitted through these n = p system points. 
The inset separately shows the contributions in eqn. (70) from the conformational 
friction corrections, cp — 2APn/a (solid lines), and from the homogeneous friction, 
ck' — 8K'n/a (chain curves). 

Figure 11: Tracer diffusion constants normalized by the Rouse result versus tracer 
molecular weight for tracer motion through gels with different pore sizes; the correla- 
tion lengths agree with the radii of gyration of tracers of the degrees of polymeriza- 
tion, N/Nf. = 10, 50, 100, 500, 1000. Tracer mass scaling exponent for a self avoiding 
walk polymer, v — 3/5, and gel fractal dimension. Dp — 5/3, are used while fixing the 
interaction strength to \g — 100 in part (a) and = 20 in (b). For a soft, polymeric 
gel, eqn. (77), the ratio of screening and entanglement lengths is fixed to ^/6 = 0.3, 
and thick solid lines are drawn. Chain curves correspond to the hard, structurally 
rigid gel case, eqn. (76). The thin line is the asymptote DN which meets 

the curves for N/N^ ~ 10^. The maximal intermediate slope, DN ~ ]\i-{4:-2u{d-Dp)) ^ 
is shown as a dotted line. In part (a), power laws, DN ~ N"^-'^ (short dashes) 
and DN ~ N~^-^ (long dashes), are compared to the soft and hard gel calculation at 
= Rg{N = 500A^e) respectively. In part (b), power laws, DN ~ N~^-^ (short dashes) 
and DN ~ iV~^-^ (long dashes), are compared to the soft and hard gel calculation at 
Cg ^ Rg{N ^ 500Ne) respectively. 
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